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A METHOD FOR DETERMINING THE TOTAL DRAG 
OF A POINTED BODY OF REVOLUTION IN SUPERSONIC FLOW 
WITH TURBULENT BOUNDARY LAYER 

By Charlie M. Jackson, Jr., and Rudeen S. Smith 
Langley Research Center 

SUMMARY 

An approximate mathematical model of the inviscid and viscid flow field about a 
body of revolution has been developed by utilizing modifications of the second-order 
shock -expansion theory and of a theoretical analysis of compressible boundary -layer 
growth with pressure gradient. The area of application of this model is intended to be 
restricted to the attached bow-shock condition in the supersonic flow region with an all- 
turbulent boundary layer. The viscid-inviscid interaction is accounted for by iteration 
between the respective theories. The solution of the mathematical model including the 
iteration and the integration of surface conditions to give wave and friction drags has 
been programed for digital computation, and the program is included herein. 

The inviscid solution is shown to compare well with the method of characteristics 
for a wide range of body shapes including local surface compressions and moderate 
amounts of boattail. The basic theory has been extended to provide an estimate of the 
flow -field pressures and shock shape. 

The viscid solution in its simplified flat -plate and conical forms is demonstrated 
to agree well with the established flat-plate and cone solutions for a turbulent boundary 
layer. The theoretical estimates of the local skin friction were in good agreement with 
experiment for a body with moderate pressure gradient; even for a condition of strong 
favorable and adverse gradients on the same body, fair agreement was indicated. 

The viscid and inviscid solutions including interactions were applied to several 
body shapes over a wide range of test conditions and the resulting total -drag estimates 
generally agreed with the reported experimental values to within 5 percent. In view of 
the difficulty involved in obtaining precise experimental data, it is believed that the pres- 
ent method is generally as accurate as the current experimental testing techniques. 



INTRODUCTION 


Currently many theoretical methods are available for predicting the viscid and 
inviscid drag forces on bodies of revolution. Most of these methods are based on 
assumptions which restrict the area of applicability. In the past two decades the effort 
expended on the development of drag prediction techniques has been primarily directed 
toward slender configurations. Engineering estimates can be obtained for the total drag 
of slender bodies of revolution in the low supersonic speed range with relatively little 
effort by applying one of the methods based on linear theory such as the supersonic area 
rule of reference 1 to obtain the wave drag and by using the turbulent -boundary -layer 
flat-plate friction analysis of reference 2 to obtain the friction drag. As the Mach num- 
ber increases or the body slenderness decreases — that is, as the hypersonic similarity 
parameter (ratio of the free -stream Mach number to the fineness ratio of the body) 
increases - the results obtained from these methods become unreliable. A need there- 
fore exists for a reliable method with which an engineer can, with relative ease, obtain 
estimates of the total drag of bodies of revolution over a wider range of conditions. 

There are, of course, well-documented theories capable of accurately predicting the 
wave drag of bodies of revolution at high values of the hypersonic similarity parameter 
(refs. 3 and 4) and theoretical turbulent -boundary -layer analyses capable of providing 
accurate estimates of the boundary -layer growth and skin -friction drag (refs. 5 and 6). 

The purpose of the present work is to modify and combine the most promising of these 
theories to obtain a mathematical model of the viscid and inviscid flow field with sufficient 
detail to account for the most important interaction effects. 

In the present paper the mathematical model is developed by using a modified form 
of the second-order shock-expansion theory of reference 3 to define the inviscid flow 
field. The turbulent -boundary-layer growth characteristics are obtained from a modifi- 
cation of the theory of reference 5, and an iteration technique is used to provide the final 
viscid and inviscid flow field from which the total drag is obtained. For convenience, 
the method is programed for the Control Data 6600 computer system and the program 
listing is included. 


SYMBOLS 

A cross-sectional area of body 

a sonic velocity 

a* sonic velocity at T* 
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C n drag coefficient, , a &-— 

a A 

max 

C p average coefficient of friction based on body surface area and free-stream 

conditions 

C p pressure coefficient, — — — 

^oo 

Cj local coefficient of friction based on free-stream conditions 

FR fineness ratio 

f friction term defined by equation (A5) 

K slope of conical surface 

l total body or flat -plate length 

M Mach number 

N velocity profile index (N = 7 in present analysis) 

n constant (n = 1 for laminar flow; n = 1/6 for turbulent flow) 

p static pressure 

q dynamic pressure 

R transformed body radius 

Rg gas constant 

R m unit Reynolds number based on free-stream conditions 

R , free-stream Reynolds number based on body length 

r local body radius 
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transformed longitudinal distance along body surface from nose 


longitudinal distance along body surface from nose 

temperature 

reference temperature 

velocity 

longitudinal coordinate (on body center line) 
radial coordinate (origin at body center line) 
transformed radial coordinate (origin at body surface) 
radial coordinate (origin at body surface) 
shock angle (see eq. (1)) 

ratio of specific heats, 1.4 for present analysis 
boundary -layer thickness 


displacement thickness, J* ^1 - + ^-jdy' 

streamline displacement due to boundary layer (see eq. (17)) 
half -angle of conical frustum body element 


entropy 


momentum thickness, 



dy' 


absolute viscosity 


kinematic viscosity, p/p 



p 


density 


Subscripts: 

base model base 

c cone 

e conditions at edge of boundary layer 

exp experimental 

i incompressible 

inv inviscid 

lam laminar 

max maximum 

T transition point 

t stagnation conditions at edge of boundary layer 

turb turbulent 

w wall 

°° free-stream conditions 

A bar over a quantity refers to flat -plate conditions. 

THEORETICAL METHOD 

The purpose of the present paper is to develop a mathematical model of the viscid 
and inviscid flow in the region of the surface of a body of revolution. This model is con- 
structed with a view to providing an accurate total -drag prediction technique for arbitrary 
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bodies of revolution in the supersonic speed range with turbulent boundary layers. In 
order to construct the model, it is necessary to provide an inviscid solution and a viscid 
solution such that the two can be iterated to provide interaction effects. In the present 
paper these solutions are provided by modifying existing theories, as described in the 
following sections. 


INVISCID SOLUTION 

The inviscid solution of the present analysis is based on the second-order shock- 
expansion theory presented in reference 3. The details of the theory are well docu- 
mented and are not discussed herein. Instead, a general description of the techniques 
used to employ the theory is given. 

In order to apply the second-order shock-expansion theory, the body is repre- 
sented by a series of cone frustums as shown schematically in figure 1. Also shown in 
figure 1 is a typical surface -pressure distribution obtained for the body by the second- 
order shock-expansion methods. The starting flow conditions at point 1 are based on the 
conical flow generated by the initial body section and may be obtained from the tables of 
reference 7. At the juncture of the initial cone and the first frustum element, the pres- 
sure change is determined by a Prandtl-Meyer expansion. The variation of the pressure 
on the frustum element from 2 to 3 in figure 1 is obtained by evaluating the first few 
terms of a series which would allow the surface pressure to approach conical pressure 
if the frustum were infinitely long. With this series, the pressure at 3 is obtained; then 
another Prandtl-Meyer expansion is used to obtain the pressure at 4; and so forth. The 
present technique is to divide the body into many frustums and base the body pressure 
distribution on the pressures at the midpoint of each element. An inherent difficulty 
with the second-order shock-expansion method is pointed out in reference 3. This dif- 
ficulty arises at higher Mach numbers when the series which describes the pressure 
variation over an element becomes divergent. When this condition occurs, the present 
technique simply assumes the pressure to be constant over the element and equal to the 
value obtained after the expansion from the previous element. 

When the second-order shock -expansion theory is applied to the boattail region of 
a body, the pressure series cannot be made to approach a conical value. In this region 
the present analysis requires the pressure series to approach a free-stream value at 
infinity, as suggested in reference 3. Figure 2 presents a comparison of the surface- 
pressure coefficient Cp determined from the present analysis (body represented by 
20 cone frustums) with that determined by the method of characteristics for typical boat- 
tail bodies (Sears-Haack bodies) of fineness ratio 4 and 10 at M^ = 2.0. The results 
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presented in figure 2 indicate excellent agreement of the second-order shock-expansion 
theory with the method of characteristics except on the aft portion where the boattail 
angle exceeds 10°. 

Since the viscid solution of the present analysis will be based on inviscid conditions 
at the edge of the boundary layer, it is important that estimates of the pressure and Mach 
number be obtained in the region of the body surface. In order to conveniently provide 
these estimates, it has been assumed that the surface pressures and Mach numbers are 
propagated out from the surface along local Mach lines. For example, in figure 1 the 
surface pressure and Mach number in the region 1-3-b-a is assumed constant and equal 
to the surface pressure and Mach number at point 2-3. The boundaries of this region, 
lines 1-a and 3-b, are assumed to be Mach lines generated by the surface Mach numbers 

at points 1 and 3, respectively. Thus the angle of the boundary 1-a is sin - l -i-+ e n ■, 

Mi u-i 

and the angle of the boundary 3-b is sin - l — p + 6 2-3' P ressures and Mach numbers 

O 

throughout the entire flow field, including the bow-shock shape, can be estimated with 
these assumptions. The shock angle can be determined for each constant -pressure ele- 
ment of the flow field from the equation 


ot = sin" 1 \| °° 9 (1) 

V 7M J 

An example of this procedure is shown in figure 1. In figure 1 the shock angle for the 
starting cone is obtained from the tables of reference 7 along with the surface Mach num- 
ber. The surface Mach number is used to obtain the direction of the Mach line from the 
surface of the body at 1. The intersection of the starting cone shock with the Mach line 
from 1 defines the shock point a, from which the next shock element is constructed. The 
shock element from a to b is constructed at the angle a 2 -3 w ^ich is computed from 
equation (1), with the pressure behind the shock defined as the surface pressure P 2 _g. 

In this manner the preliminary estimates of the shock shape and flow-field pressures 
are determined for the entire flow field. 

In order to evaluate the estimates of the flow -field pressures and shock shape, 
results of the present method and the method of characteristics are compared in figure 3. 
In figure 3 the pressure coefficients are shown for Sears-Haack bodies of fineness ratios 
4 and 10 at a free-stream Mach number of 2.0. Good agreement is indicated between the 
present method and the method of characteristics in the region of the body surface, with 
the exception of the aft portions of the body where high boattail angles occur. These 
results also indicate that the present method can be used to predict the shock shape and 
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forebody flow-field pressures to sufficient accuracy for preliminary interference 
analysis. 

The second-order shock-expansion methods of reference 3 can be applied to bodies 
with local compressions, provided the compressions are small enough for the shock to be 
theoretically attached to the compression corner. The details of the calculating proce- 
dure for compression elements on the body surface are given in reference 3. The pro- 
cedure essentially involves the use of the two-dimensional oblique shock relations at the 
start of the element and the evaluation of a series to provide the pressure variation over 
the element in the same manner as described for an expansion element. Reference 3 
did not illustrate the use of the second-order shock- expansion theory on a body with local 
compression, so in order to evaluate the theory under these conditions, the pressure dis- 
tribution on an inlet spike has been calculated by this method and compared with that cal- 
culated by the method of characteristics. This comparison is shown in figure 4 for a 
free-stream Mach number of 3.0. Figure 4 indicates good agreement between the theories 
for both the compression and the following expansion. For the present method the inlet 
spike of figure 4 was represented by 19 elemental cone frustums. 

VISCID SOLUTION 


The viscid solution of the present analysis consists of the use of the integral 
momentum equation in conjunction with a mathematical model of the turbulent boundary- 
layer growth. The present theory for determining the turbulent boundary -layer growth 
is based on a modification of the analysis of reference 5. The mathematical expression 
for the variation of the momentum thickness as developed in reference 5 is (in the nota- 
tion of the present paper): 


e = 


(l + 0.2M e 2 ) 3 

0.0076| — 
\ a t, 

r 

M e 3 r 


-5/42 


^1 + 0.2Mg ^ 

1 



M e 10 /3 r V6 
(l + 0.2M e 2 ) 4 


-I 6/7 


ds 


( 2 ) 


Equation (2) is developed in reference 5 by extending the incompressible equation of ref- 
erence 8 to the compressible condition by a transformation based on the work of refer- 
ence 9. The analysis of reference 5 was slightly modified as indicated in appendix A, 
and equation (2) was rederived for the purposes of the present analysis to give 

16/7 


e = 


(l + 0.2M e 2 ) 3 

r . . 1/6 

v t\ 

0.0076 -i 

V t) 

p M e 10 / 3 r 7 / 6 

M e 3 r 

( 2) -11 / 42 1 
(l + 0.2Mg j J 

0 (l + 0.2M e 2 ) 4 


ds 


(3) 
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A comparison of equation (3) with equation (2) indicates the net effect of the modifications 
of the present analysis to be the alteration of an exponent of a term in the denominator. 
This change has a tendency to increase the magnitude of the estimated value of the 
momentum thickness for specific inviscid conditions. 

The effects of the modifications of the present analysis on the accuracy of the 
momentum -thickness expression of reference 5 are indicated in figure 5. Figure 5 
shows the ratio of the momentum thickness calculated from equations (2) and (3) to the 
momentum thickness obtained from velocity profile measurements reported in refer- 
ences 10 and 11. In general, better agreement with the available data is indicated for 
the modified equation of the present analysis (eq. (3)) than for equation (2). Therefore, 
equation (3) is used herein. 

A mathematical model of the turbulent boundary -layer growth on a body of revolu- 
tion can be obtained on the basis of equation (3) with the additional assumption of an ana- 
lytical form for the velocity profile. In the present analysis the well-known power -law 
profile is assumed to represent the velocity profile shape. Expressed in the present 
notation, the power -law profile is 


u 




( 4 ) 


where the value of N is taken to be 7 for all conditions. The analysis of reference 12 
incorporates equation (4) to provide expressions for the boundary -layer thickness and 
displacement thickness as functions of the momentum thickness on bodies of revolution. 
The expressions of reference 12, in the notation of the present paper, are 
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and 



In the present analysis equations (5) and (6) are used in conjunction with the momentum 
thickness obtained from equation (3). By an iteration process the boundary -layer param- 
eters 5, 6*, and 6 are obtained as functions of the body radius and inviscid flow con- 

ditions at the edge of the boundary layer. In order to start the iteration the initial values 
of 9/6 are taken to be those presented in the tables of reference 13 for two-dimensional 
boundary layers. 


The local skin friction can be obtained from the variation of the boundary -layer 
parameters provided by the expressions of equations (3), (5), and (6) and the momentum 
equation. The form of the von Karman integral momentum equation appropriate to bodies 
of revolution with pressure gradient is 


^f 0* _ d0 + q 1 dr | 
2 q e ds r ds 


2 - M e 2 + — 
e 9 

M e (l + 0.2M e 2> 


dMe 

ds 


(7) 


Before considering the details of the solution of the boundary -layer equations for an 
arbitrary body with pressure gradient, it would be interesting to examine the expressions 
of the present analysis for flat-plate and conical flows for which closed-form solutions 
are available. Equations (3) and (7) can easily be reduced to flat -plate conditions by 
setting M e = and q g = q^ and by letting r — °°. The following results for the 
flat -plate momentum thickness and local skin-friction variations are obtained: 

j 0.0152s 6 / 7 ft 1 , j 

(l + 0.2Mj) 0 ' 41 V/ 7 V T t +216 /' 


( 8 ) 


( 9 ) 



0.0261 

2\0.418, 


/ 2W- 418 / > 1/7 

(1 + 0.2MJ) (R„s ) L/ ‘ 


Too + 216 

T t + 216 



In these equations the Sutherland law of viscosity has been used to relate the viscosity at 
stagnation conditions to that at free -stream conditions. The equation for average flat- 
plate friction can be obtained from equation (9) by integrating Cf from 0 to l to yield 



0.0306 


,1 + 0.2M C 


,0.418 
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1/7 


/T 0O + 216\ 1/7 
\T t +216 ) 


( 10 ) 


The basic equation for this analysis (eq. (3)) was developed by using the incompres- 
sible average skin-friction law of Falkner and a compressibility correction based on an 
approximation to the method of reference 2 (see appendix A). This compressibility cor- 
rection can be obtained by dividing equation (10) by the Falkner law, with the following 
result: 



TOO+ 216\ 1/7 

T t + 216 / 


,0.418 


1 + 0.2M„ 


(ID 


A comparison of the results of equation (11) with the results of the Sommer and 
Short T' method of reference 2 is presented in figure 6. For the method of reference 2, 
the Karman-Schoenherr incompressible friction law and a Reynolds number of 10 X 10® 
were used to obtain the compressibility factor. Figure 6 indicates that equation (11) 
provides a simple approximation to the Sommer and Short T' method of reference 2 with 
a fair degree of accuracy for the conditions of interest in the present analysis. 

The results for turbulent flat -plate skin friction obtained by the Sommer and 
Short T* method of reference 2 are compared with the results of equation (10) in fig- 
ure 7. The comparison indicates fair agreement between the methods for the conditions 
shown. However, skin friction obtained by the present method appears to vary much less 
with Reynolds number than does skin friction obtained by the Sommer and Short T' method. 
This difference is due to the basic difference between the Falkner law used in the present 
analysis and the Karman-Schoenherr law used in the Sommer and Short T' method. 


The basic equations of the present viscid analysis (eqs. (3) and (7)) can be reduced 


for the condition of conical flow by setting 


%_ = M e p e 
q oO M 2 Poc 


■, r = Ks, and 


dM ( 

ds 


= 0. The 
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resulting expressions for the momentum thickness and local skin -friction coefficient on 
a cone are 


0 c " 


0.00787 exp^ 

Av/lRg) 

M„1/7 s 6/7 ( 

(T^+216 

M e Roo(l + 0.2MJ 2 

^1.5~ 

1/7 

(l + 0.2M e 2 ) 

^0.204' 

V T t + 216 


,1 n 


( 12 ) 


and 


V - 0.0292 


M\ 13/7 (i + 0.2M oo 2 ) 3 ‘ 286 /T oo+ 216 ! \ 1/? 

(l + 0.2M e 2 ) 3 '^ 04 T ' + 216 R °° ' 


(13) 


The local coefficient of equation (13) can be integrated over the surface of the cone 
and the resulting average friction coefficient can be written as follows: 


/ 6 Ar?\/ M e\ 3/7 ( 1 + °’ 2M «> 2 ) / T oo +216 1 \ 7 

C„ „ = 0.0315 exp (- — j ( — | 

F ’ c \ 7R gA M »/ / 9 v3.704\T t + 216 R* J 

(l + 0.2Mp 2 ) ' 1 ,l/ 


(14) 


A comparison of these expressions with flat -plate expressions can be made by 
dividing equations (13) and (14) by the flat -plate expressions of equations (9) and (10), 
respectively. The results of this operation are 



1.120 exp 


6 A n \/ M e \ 13 / Y 1 *°- 2M «' 2 

7 Rg J\Moo/ yi + 0.2M e 2 


(15) 


and 


. ... .13/7/ 

2 ! £= 1.035 e X p('4^Y— ) 

C F V 7R gA M »/ \ 


1 + 0.2M k 
1 + 0.2M e 2 


,3.704 


(16) 


Ah examination of equations (15) and (16) indicates that, as the cones become slender, 
M g — and Ari/Rg — 0 which reduces c^ c jc ^ and Cp c ^Cp to 1.120 and 1.035, 
respectively. Reference 10 gives representative values for these ratios to be 1.194 and 
1.04, respectively, on the basis of the work of Van Driest and others. 


12 



VISCID -INVISCID INTERACTION 


In the present analysis it is assumed that the viscid -inviscid interaction is suffi- 
ciently weak for the flows to be calculated separately, as outlined in the previous sections, 
and the interaction effects determined by iteration. The details of this technique can best 
be understood from an example problem. As a typical example, the calculation of the 
viscid and inviscid flow over an Adam’s body with fineness ratio 10 and with a ratio of 
maximum area to base area of 1.88 is considered. Free-stream flow conditions are 
taken to be = 4.0 and ^ = 9.0 x 10®. The variations of the important param- 
eters of the first two iterations for the inviscid and viscid solutions are presented in fig- 
ures 8(a) and 8(b), respectively. Initially, the inviscid analysis of the second-order 
shock-expansion method is applied to the body shape. In this example the body was rep- 
resented by 70 elemental cone frustums. The resulting variations of Cp and M e are 
shown in figure 8(a). With this initial inviscid solution assumed to provide conditions at 
the edge of the boundary layer, numerical solutions of the viscid equations can be pro- 
vided and the results yield the boundary -layer growth and skin-friction coefficients 
shown in figure 8(b) for the first iteration. At this point a first-order mathematical 
model is available to describe the inviscid surface conditions and the boundary -layer 
growth. 

The accepted method of evaluating the weak viscid-inviscid interaction is to calcu- 
late the inviscid solution for an effective body shape which is obtained by including the 
displacement thickness of the boundary layer. The boundary -layer growth is then reeval- 
uated according to this effective -body inviscid solution. In the present analysis, the effec- 
tive body shape is obtained from the displacement thickness as suggested in reference 14. 
The surface streamline displacement due to the presence of a boundary layer on a body 
of revolution in reference 14 is related to the displacement thickness of the present 
analysis by the following equation: 



The effective body shape obtained for the second iteration of the example problem 
is shown in figure 8(a) along with the variations of Cp and M e for this body shape. 
The effects of the viscid solution on the inviscid solution are indicated by the differences 
between the first and second iterations shown in figure 8(a). In order to obtain the effect 
of the change in the inviscid solution on the viscid solution, the results of figure 8(a) are 
used to reevaluate the boundary -layer growth and skin friction. The results of the 
reevaluation are shown as the second iteration in figure 8(b). 
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If the boundary layer is considered to be of finite thickness, then the inviscid con- 
ditions on the surface of the effective body do not represent those at the edge of the 
boundary layer. In the present analysis it is assumed that the Mach number on the sur- 
face of the effective body is propagated out to the edge of the boundary layer along a Mach 
line. 

This assumption is represented schematically in figure 9. The geometry of fig- 
ure 9 indicates that the Mach number distribution at the edge of the boundary layer can 
be approximated by the variation of the Mach number calculated for the surface of the 
effective body and transformed by the following relation: 

Ax * (6 - 6**))jM im 2 - 1 (18) 

In equation (18) the value of 6 is obtained from the previous iteration. 

The Mach number distribution of figure 8(a) (second iteration) was transformed by 
equation (18) and used in the solution of the viscid equations to provide the second- 
iteration solution of the boundary -layer growth and the attendant skin friction. The com- 
parison of the first- and second-iteration boundary-layer growth characteristics shown 
in figure 8(b) indicates little change. Further iteration did not yield significant changes 
in the solutions, and the second-iteration values of figure 8 are considered to be repre- 
sentative of the viscid -inviscid mathematical model of the flow including interations in 
the region of the body surface. The wave and friction components of the total drag can 
then be obtained by integration of the pressure and local friction coefficients over the 
body surface. 

It has been mentioned previously that the solutions of the boundary -layer equations 
of the present analysis were numerical in nature and required the use of high-speed 
digital computers. Thus, the entire process described in the solution of the example 
problem has been programed for high-speed digital computation. A detailed description 
of the resulting computer program is presented in appendix B along with the necessary 
operating details. The program is designed to make as many iterations between the 
inviscid and viscid solutions as the user desires. For the remainder of the present 
analysis, all theoretical solutions, unless otherwise noted, will be the results of the sec- 
ond iteration obtained from the computer program presented in appendix B. 

One difficulty which was encountered in the use of the present method involves the 
selection of the proper starting cone for the inviscid solution. This difficulty occurs 
when the shape of the body is such that a wide range of starting cones can be made tan- 
gent to the body surface in the region of the nose (i.e., a body with a rapid change of slope 
in the nose region). For bodies of this type it is recommended that the starting cone be 
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chosen such that the boundary layer is within the surface stream tube generated by the 
flow field of the starting cone. If this condition is satisfied, then for the expanding flow 
the entropy at the edge of the boundary layer remains constant and equal to the entropy 
of the flow field of the starting cone (an assumption of the second-order shock-expansion 
theory). An approximate criterion for the condition that the boundary layer be within the 
surface stream tube can be expressed as 

y a 2 > [(6 + r) 2 - r 2 ] (19) 

where y a is the radial coordinate of the conical shock point a in figure 1. 

A criterion for the selection of the size of the body elements would be desirable. 
Unfortunately, no criterion is available to cover the wide range of body shapes and Mach 
numbers available to the user. It is suggested that the body be described by no less 
than 20 elements including the starting cone and, if the body shape is irregular (i.e., 
adverse and favorable gradients in combination), that this number of elements be 
increased until the surface pressures and skin friction show no significant change. 

EVALUATION OF THE THEORETICAL METHOD 

In this section the mathematical model which was developed in the previous section 
will be applied to several body shapes to describe the viscid and inviscid flow in the 
region of the body surface. The resulting viscid and inviscid drag estimates are com- 
pared with experimentally determined values where available. The viscid-inviscid inter- 
action effects are accounted for herein by using the second -iteration values of the local 
surface conditions to determine the drag. 

The body shape used in the first example is that of the NACA RM-10 missile body 
described in reference 10. The estimates of the momentum thickness and local skin- 
friction coefficients are shown in figure 10 along with the body shape. Figure 10 also 
includes some experimentally determined values of momentum thickness and local fric- 
tion coefficients obtained from reference 10, and a comparison of the experiment and 
theory indicates good agreement of the present estimates. 

Figure 11 compares the friction-drag and total -drag estimates of the present 
method with experimental values for the RM-10 missile body at several Mach numbers 
and Reynolds numbers. It should be noted that the experimental total-drag data of ref- 
erence 10 were obtained by using the method of characteristics to provide the inviscid 
drag contribution and adding this directly to the friction drag which was obtained from 
velocity profile measurements. The experimental data from reference 15, however, 
were obtained by force and pressure measurement and are probably more accurate. A 
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comparison of the theory and experiment shown in figure 11 indicates that the present 
theoretical method is in good agreement with the experimental data from both references. 

In order to examine boundary -layer behavior in the presence of strong favorable 
and adverse pressure gradients, reference 16 has provided experimental data on a low- 
fineness-ratio waisted body at several Mach numbers up to =2.8. The experimental 
values of the surface-pressure coefficients and local skin-friction coefficients are pre- 
sented in figure 12 along with the body shape. The theoretical estimates of the surface 
pressures and skin friction are also shown in figure 12. It should be noted that these 
theoretical estimates were obtained with five iterations between the viscid and inviscid 
flow fields. The estimated surface pressures show good agreement with the experimental 
values of reference 16 except for x/l > 0.7, and the agreement of the skin-friction esti- 
mates and experimental values is fair. One obvious difficulty is the presence of a lim- 
ited laminar run on the nose of the body (indicated by the lower values of the experimen- 
tal data in the region of the nose). Although there is no inherent reason that the present 
technique cannot be extended to include a laminar rim (the analysis of appendix A devel- 
ops the growth of the momentum thickness including an initial laminar run), the present 
analysis does not account for this condition. 

The surface pressures and local skin-friction data presented in reference 16 were 
integrated to give the friction and wave drag of the waisted body of revolution at the 
supersonic Mach numbers. The friction-drag and total -drag values are presented in 
figure 13 for a Reynolds number of 10 x 10 6 . The theoretical estimates provided by the 
present method are indicated to be in good agreement with the experimental total drag. 

The method of the present analysis was used to obtain theoretical estimates of the 
total drag of several body shapes over a wide range of flow conditions. The ratio of the 
theoretical drag coefficient predicted by the present method to the experimental drag 
coefficient is shown in figure 14 as a function of the hypersonic similarity parameter 
(see solid symbols). Also shown is the ratio of the theoretical to experimental drag 
coefficient, with the theoretical values obtained by using the supersonic area rule of ref- 
erence 1 along with the flat -plate skin-friction method of reference 2 (see open symbols). 
The experimental values used were taken from references 10, 15, 16, and 17 and from 
present data obtained in the Langley Unitary Plan wind tunnel. It is significant to note 
that the supersonic -area-rule results are not reliable for conditions at M^/FR > 0.2 
whereas the present method gives reliable results at all conditions shown. In fact, all 
the cases considered in figure 14 indicate the present method to be within 10 percent of 
the reported experimental drag values, with all except four points being within 5 percent. 

In view of the difficulties associated with obtaining accurate experimental drag data (i.e., 
attaining fully turbulent flow, measuring base pressure, etc.), it is believed that the pres- 
ent method is generally as accurate as the current experimental testing techniques. 
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CONCLUDING REMARKS 


An approximate mathematical model of the inviscid and viscid flow field about a 
body of revolution has been developed by utilizing modifications of the second-order 
shock-expansion theory and of a theoretical analysis of compressible boundary -layer 
growth with pressure gradient. The area of application of this model is intended to be 
restricted to the attached bow-shock condition in the supersonic flow regime with an all- 
turbulent boundary layer. The viscid -inviscid interaction is accounted for by iteration 
between the respective theories. The solution of the mathematical model including the 
iteration and the integration of surface conditions to give wave and friction drags has 
been programed for digital computation, and the program is included. 

The inviscid solution is shown to compare well with the method of characteristics 
for a wide range of body shapes including local surface compressions and moderate 
amounts of boattail. The basic theory has been extended to provide an estimate of the 
flow -field pressures and shock shape. 

The viscid solution in its simplified flat -plate and conical forms is demonstrated to 
agree well with the established flat -plate and cone solutions for a turbulent boundary 
layer. The theoretical estimates of the local skin friction were in good agreement with 
experimental values for a body with moderate pressure gradient; even for a condition of 
strong favorable and adverse gradient on the same body, fair agreement was indicated. 

The viscid and inviscid solutions including interactions were applied to several 
body shapes over a wide range of test conditions, and the resulting total drag estimates 
generally agreed with the reported experimental values to within 5 percent. In view of 
the difficulty involved in obtaining precise experimental data, it is believed that the pres- 
ent method is generally as accurate as the current experimental testing techniques. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 6, 1968, 

722-01-00-02-23. 


APPENDIX A 


MODIFICATION OF EXISTING EXPRESSION 
FOR MOMENTUM THICKNESS 


The analysis of reference 5 develops an equation which describes the momentum 
thickness on a body of revolution with pressure gradient in compressible flow. The 
approach taken in reference 5 is to transform the incompressible equation of reference 8 
by means of the transformations described in reference 9. The present development of 
the momentum-thickness expression follows the approach of reference 5 with these two 
modifications: 


(1) The reference temperature equation of reference 2 is used. 

(2) The viscosity is assumed to vary as the square root of the temperature. 


The starting point for the present development is taken to be the incompressible 
momentum-thickness equation of reference 8 which is expressed as (in the notation of 
the present paper) 


0i = 




u, 


, i 5 R 2 dS 


-17/12 


+ f turb^t 1//6 


S 

S T 


,6/7 


u, 


^10/3 r 7/6 


e,i 


(Al) 


The transformation of this equation to the compressible form is accomplished according 
to reference 5 by the following equations: 


-\ 


3y-l 

/ a p \ y- 1 
dS = ( —I ds 

\ a t / 

dY ' = $£ d y’ 

H 

Ui = — u 


u e,i = a; u e 


(A2) 


R = r J 

In the process of this transformation, reference 5 and the present analysis make the 
assumptions that there is no heat transfer, that the pressure is constant through the 
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APPENDIX A - Continued 


boundary layer in the direction normal to the body surface, that the Prandtl number is 
unity, and that the flow external to the boundary layer is isentropic. 


The left-hand side of equation (Al) can be conveniently transformed by equa- 
tions (A2) if, for purposes of this transformation only, the boundary layer is assumed 
to be small with respect to the body radius so that the momentum thickness can be 
expressed as 



(A3) 


With the assumption indicated in equation (A3) and the transformations of equations (A2), 
the incompressible momentum thickness is related to the compressible value by 

y + 1 
y - 1 

9 (A4) 



Most of the terms on the right-hand side of equation (Al) can be directly trans- 
formed by equations (A2). However, the transformation of the f terms is complicated 
by the fact that according to reference 5 they are related to the average flat-plate fric- 
tion coefficients by 


f = 


/ •, \ n/ — \ 1+n 

u e,i^ | CF,ij 


(A 5) 


The value of n is 1 for laminar conditions and 1/6 for turbulent conditions. 


In order to transform the f terms, it is first necessary to establish the trans- 
formation for the average flat-plate friction. Reference 5 shows that if the Blasius equa- 
tion is used for the laminar average friction relation, the transformed value of fj ain 
does not differ from the incompressible value which is determined to be 0.441. The value 
of ft ur b on the other hand is shown to be 0.0076 for the incompressible case when the 
Falkner relation is used for CF,b and the transformation to the compressible case 
involves the reference temperature concept. The present analysis follows the trans- 
formation of reference 5 with a modified reference temperature relation and a modified 
viscosity-temperature relation. The reference temperature equation of the present 
analysis is that of reference 2 for adiabatic wall conditions and is given as 


rp* O 

4r-= 1 + 0.1142Mg = 
A e 



The viscosity-temperature relation of the present analysis is 


Me 



a e 

at 


(A6) 


(A7) 
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APPENDIX A - Concluded 

The average flat-plate incompressible skin-friction coefficient obtained from the 
Falkner equation can be transformed by equations (A2), (A6), and (A7) to provide the 
expression for the compressible skin-friction coefficient as follows: 


aeWVM 11 / 7 0.0306 


P t u i z i\ 1/7 


F,turb l at J la* 


The curve fit of figure 15 shows that equation (A8) can be further simplified without sig- 

/ /a \5A 

nificant loss of accuracy in the 0 s M e = 5.0 range by replacing a e/ /a* with ( — ) 


With this simplification and equations (A2) and (A5), the results of equation (A8) can be 
combined to provide the desired transformed coefficient f^^. The value of f^^ can 
be simply expressed as 


/at \ 11/21 

f^b = 0.0076U 


With the results of equations (A2), (A4), and (A9), the incompressible momentum- 
thickness expression of equation (Al) can be transformed to the compressible form with 
the following result: 


( 1+2 TT M e 2 r n f S T M. 5 r 2 ds 
e = ± w o.44i § 


M e 3 r 


0.0076^) 1/6 


1 + “ — — — ‘ M, 


1 „ o\4 


M l°/ 3 r' 7 / 6 ds 


i + ^Me 2 )" 11 ^ Js T (l + ^M e 2 


(A10) 


For the condition of an all-turbulent boundary layer, equation (A10) becomes 


**r o- 0016 ^ 176 


/ 9\3 0.0076 — 

_ jl + 0.2M e 2 j \*t 


10/3 r 7/6 ds 


M e** r (l + 0.2M e 2 ) Jo (l + 0.2M e 2 ) 


(All) 
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COMPUTER PROGRAM TO OBTAIN VISCID AND INVTSCID DRAG 
CONTRIBUTIONS FOR A BODY OF REVOLUTION 

A process is discussed in the text for obtaining the inviscid solution of the flow field 
and the solution of a mathematical model of the boundary-layer growth on a body of rev- 
olution. This process includes the numerical solution of the viscid and inviscid equations 
as well as an iteration between the two solutions to provide interaction effects. The 
numerical solution and iteration have been programed for high-speed digital computation, 
and the program includes an integration of the local surface pressures and skin friction 
to provide the respective drag contributions. The purpose of this appendix is to provide 
a description of the necessary input and available output as well as a FORTRAN listing of 
the source program with an example input case and the resulting output listing. 

Description of Program 

The program reads the body geometry in terms of x,y coordinates which are used 
to represent the elemental conical frustums of the body in the numerical process. The 
starting cone is generated by extending the slope of the first elemental frustum forward 
to provide the full cone. The inviscid solution is then calculated for the midpoints of the 
frustum elements. Numerical derivatives are taken from linear interpolation of the 
inviscid quantities, and the viscid equations are solved. 

The boundary -layer growth is then used to provide an effective body shape for the 
second iteration. The program continues iterating between the viscid and inviscid 
solutions for a specified number of iterations. The inviscid flow fields and the resulting 
boundary -layer growth are output for each iteration along with the wave and friction 
drags. For convenience the flat-plate friction drag is output along with the values for 
the body of revolution. 


Program Listing 

The FORTRAN listing of the source program used at the Langley Research Center 
on the Control Data 6600 computer system is presented as follows: 
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PROGRAM V8GDY ( INPUT , OUTPUT , TAP E5 = [ NPUT , TAPE6 = OUTPUT I 
DIMENSION WHAT(ld) , X(200», Y(2D0), XANSI200), PANS(200), AM ANS ( 20 
10), CP ( 20U ) , XS( 200 ) , YS ( 200 ) , DEL(200), ANGLE ( 200 ) , PC(200), A<4) 
2, TEMPI8), ROOT! 3) , DL SOR ( 1 3 , 1 1 ) , DELT AT ( 13 , 1 1 ) , PCT(13,11), AM1TI 
313,11), THET1TI 13,11 ) , MTAB(ll), CDWS(200), XT(200), CFTI200), CFT 
421200), MT2I200), DLSTAR(200), YS AVE ( 200 ) , THODLT(200), V(200), DE 
5LTT ( 200 ) , XT T( 200 ) , XL(200), XOLI203 ), YOLI200), THE T AT ( 200 ) 

COMMON XANS, AMANS , MT2 , DLST AR , N , XT , S, CFT ,C FT2 , V , DELSOR , AM,RN, THODLT 
1 ,X, I ,OELTT,Y ,0,YSAVE,PANS,DEL,T,CDWS,XL, THETAT.XTT, ANGLE, CP, XOL, YO 
2L,YMAX, ITER, ICUDfc, STI20C ) 

COMPLEX TEMP, ROOT 
REAL MIDDLE ,MT2,MTAB 

NAMELIST /HUM/ DEL T AT , AM , PC T , AM1T , THET 1 T , N , X , Y, MT 2 , DL SOR , RN , THODLT 
1 , 0 , M TAB , I TER , TO 

1 READ (5,29) WHAT 
READ ( 5 , NUM ) 

I CODE = l 

T=T0/( 1 , +AM**2/5. i 

DO 2 1=1, N 

XOL ( I) = X(I)/X(N) 

DEL T T ( I )=0. 

DL ST AR ( I ) =0 . 

YOL ( I ) = Y(I )/X(N) 

2 YSAVE ( I ) =Y ( I ) 

YMAX=YS AVE ( 1 ) 

DO 3 1=2, N 

3 YMAX=AMAX1(YMAX,YSAVE( I) ) 

ST( 1 )=SQRT( XI 1 ) **2 + Y( 1 )**2 ) 

DO 4 1=2, N 

4 ST( I )=SORT ((Y(I)-Y(I-l) ) **2 +■ ( X ( I ) -X ( 1-1 ) )**2 )+ST( 1-1 ) 

S=3. 1416*YMAX**2/144. 

5 WRITE (6,33) WHAT 
WRITE (6,35) ICODE , AM, RN 
WRITE (6,37) 

DEL ( 1 1=57.295779* A TAN((Y(2)-Y(1))/(X(2)-X(1>)) 

DO 6 1=2, N 

6 DEL ( I ( = 57.29577 9* ATAN( ( Y ( I ) -Y ( l - 1 > ) / ( X ( I ) -X( I -1 ) ) ) 

WRITE (6,38) (X( I),Y(I ) , DEL ( I) , 1 = 1 , N ) 

DEL ( 1 )=ATAN( (Y(2)-Y(1))/(X(2)-X(1))> 

CALL BILUP ( DELT AT , MT AR , PCT , AM IT , THE T 1 T , DL SOR , 13,11,DEL(1),AM,PC(1 
1 ) , A M A N S ( 1) ,THETA1, DELSOR) 

THET A = THET A1 *5 7.29 577 

WRITE (6,36) AMANS ( 1 ), THETA , DELSOR 

XANS(l)=X(l)/2. 

PANS ( 1 ) =PC ( 1 ) 

CP(1)=(PANS(1)-1.)/(AM**2*.7) 

AMU= AT AN ( ( 1./ AMANS ( 1 ) ) / S ORT ( 1 . - ( 1 . / AM ANS ( 1 ) ) * *2 ) )+DEL( 1 ) 

YS(1) = (Y(1)-(Y( 1 ) /TAN (DEL ( 1 I ) >*SIN( AMU ) /COS( AMU ) I /( 1. -S I N( AMU ) /COS 
1( AMU)/SIN( THET A1 )*CGS( THET A1 ) ) 

XS(1)=YS(1)/SIN(THETA1 )*C0S(THETA1 ) - Y ( 1 ) / TAN ( DE L ( 1 H * X ( 1 ) 

CDWS ( 1) =CP< 1 )*<YSAVE( 1 ) **2 ) *3 . 1416 
AMA= AMANS ( 1 ) 

ANGLE ( 1 ) = (0EL(1 ) +A S IN ( 1 . / AM A ) ) *57 ,29577 
DPDS A=0 
PA=PC ( 1 ) 

WRITE ( 6,39) 

WRITE (6,40) X ANS( 1) , PANS (1 ) , A MANS ( l),CP(l),XS(l),YS(l), ANGLE ( 1) ,C 
1 DWS ( 1 ) 
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DG 25 I =2 » N 

OtH I ) = ATAN ( CY( 1 i-Y< 1-1 I ) /( X( I l-X ( 1-1 ) > ) 

IF ( DELI 1) .LT.O. ) GO TO 7 
DELT A=DEL( l ) 

CALL BILUP (DELTAT , MT AB , PCT , PCI , PCT , PCT , 1 3 , 1 1 , DEL (I) , A M, PC ( I ) , PC ( I 
1 ) tPC t I ) »PC( I > ) 

GO TO 8 
PCI I ) = 1 . 

ANUA=2.4495*ATANt . 40 82 5*SQRT ( AMA**2- 1 . ) ) -ATAN I SORT! AMA**2-1 . > ) 
ANUB=ANUA+DEL< I-1I-DELI I* 

IF ( ANUB.GE.ANUA) GO TO 14 

IF I ABSIDELl I-l)-DEL(I) ).LT».C0301) GO TO 14 

CALCULATION FOR COMPRESSION ANO SHOCK FORMATION 

ANG= DEL I I I -DEL I 1-1 > 

A ( 1 )=1 .0 

A( 2 )=-( AMA**2+2. )/ AMA*#2-1.4*S1N(ANG)**2 

At 3 ) =< 2 . *AMA**2+1 . ) / AM A**4+ I 2 . 4**2 / 4 . + . M AM A**2 I * SI N I ANG ) **2 
A ( 4 ) =-COS ( ANG ) **2/ AM A* *4 
CALL FALG < A , 3 , G , R 00 T , TEMP , I ER R > 

AR1=REAL (ROUT { 1 ) > 

AI1=AIMAG( ROOT 1 1 ) ) 

AR2 = REAL I ROUT ( 2 I ) 

AI2 = AIMAG( ROOT (2 I ) 

AR3 = REAL (ROOT (3) 1 
A I 3 = A I MAG ( ROUT (3 ) ) 

IF ( IERR) 9,10,9 
WRITE (6,30) HOOT 
GO TO 27 

TEST TO FIND IMAGINARY ROOTS 

IF I AB SIAIl). LT. . 30 C'U L. AND. ABSIAI2). LT.. 00001. AND. AB SIAI3I.LT. ,000 
1^1) GO TO 11 

WRITE (6,31) A(2) , A(3> , A ( 4 » ,RDOT 
GO TO 27 

CHECK FOR MIDDLE ROOT 

11 B IG= AMAX1 ( AK1 , AR2, AR3) 

SMALL = AM IN1( AR1,AR2,AR3) 

DO 12 J=l,3 

IF ( REAL(ROOT (J) ). LT.B IG. AND. REAL(ROOT (J )) .GT .SMALL) GOTO 13 

12 CONTINUE 

WRITE (6,32) ROOT 
IF ( I .NE .2) GO TO 22 
PB=P A 
AMB= AMA 
GO TO 22 

13 MIDDLE 3 SORT ( RE AL( ROOT ( J ) ) ) 

SIGMA=AS IN(MIODLE) 

SINSO=SIN(S I GMA ) ** 2 

AMB= SORT ((36.*AMA**4*SINSQ-5.*( AMA*< t 2*S I NSQ- 1 . ) * ( 7 . *AM A*#2* S I NSQ +5 
1 .))/(( 7 . *AMA**2*SI NSQ-1. )*( AMA**2*SINSQ+5. ) ) ) 
PB=(7.*AMA**2*SINSQ-1. I/6.*PA 
AMUB = SIN(1./ AM B ) 

XX=( 2.4*TAN(DEL( I ) -DEL ( 1-1 * ) *COS ( S I GM A-DEL ( I - 1 ) ) -S IN( S I GMA-DEL ( 1-1 
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1 ) ) ) *AMA**2*SIN(SI3MA-DEL( I- 1 » )**2*S1N( S I GM A-DEL ( I - 1 ) 1 

YY= 1 . -t- < 1 .-2.*SIN(SI GMA— DEL ( 1-1 ) > **2 +2 . *TAN ( DEL ( I) -DEL (1-1) )*SIN( SI 
1GMA-DEL ( f-1 f ) *COS( SIGMA-OEL ( 1-1 )) ) *AMA**2*SIN { SIGMA-DEL( 1-1 ) )**2 
F=4./2.4*( l. + . 4/2. *AMA**2 )*S IN (SIGMA-DEL ( 1-1 ) )*XX/YY 
BA=1 . 4*P A*AMA**2/ ( 2.*( AMA**2-1 . ) ) 

BB=1 .4*PB*AMB**2/( 2.*( AMB**2-1 . ) ) 

XY=2 .*BB/( Y( I )— 1. I *( SIN ( SIGMA- DEL ( 1-1 ) ) *S IN (DEL (I -11) /SI N( S I GMA-DE 
1 L ( I) ) — S I N ( DEL ( I) )) +DPDS A* (BB*S INI SIGMA-DEL ( 1-1) )/( BA*SIN( SIGMA-DEL 

2 ( I ) ) )■*•( PB/PA-F)*COS( SIGMA-DEL ( 1-1 ) ) *TAN( AMUB) / S I N ( S I GM A-DEL ( I) )) 
DPDS=XY / ( 1 • + TAM AMUB ) / T AN ( S IGM A-DEL ( I) ) ) 

GO TO 21 

14 DELM = . 0 1 
K = 1 

AMB= AMA 

15 AMB= AMB + DELM 

AN U= 2. 44 95* A TAM . 40 825* SORT ( AM B** 2-1 .) ) -AT AN ( SORT (AMB**2-1. ) ) 

IF (ANU-ANJB) 15,15,16 

16 AMB=AMB-DfcLM 

GO TO ( 17, 18,19,20 ) , K 

17 DELM=.00 1 
K = 2 

GO TO 15 

18 DELM=.0001 
K = 3 

GO TO 15 

19 DELM=.D0001 
K = 4 

GO TO 15 

20 PB = P A*( (l.+AMA**2/5.)/(l.+AMB**2/5.))**3.5 
BA= 1 .4*PA*AMA**2/( 2.*( AMA** 2-1 . ) ) 

AOA STA= ( (125. /216. )/AMA)*(l , + AMA**2/5. ) **3 
BB= 1 . 4*P B*AMB**2/( 2.*( AMB**2~1 . ) ) 

AOA STB= ( ( 125./ 216. ) / AMB ) * < 1 . + AMB**2/ 5 . ) **3 

DPDS = ( BB/Y( 1-1 ))*( AOASTA/AOASTB*S IN ( DEL ( 1-1 > )-S IN (DEL ( I )) ) + BB/BA*A 
10AST A/AOASTB*DPDSA 

21 ETA=DPD S*( X ( 1 ) -X ( I -1 ) ) / ( ( PC ( I ) — P8 ) *COS ( DEL ( I ) ) ) 

IF (ETA) 22,23,23 

22 WRITE (6,34) 

PANS ( I ) = PB 
AMANS( I ) =AMB 
AMA= AMB 

P A = P B 

X AN S ( I ) - ( X ( I)+X( 1-1) )/2. 

DPDSA=0. 

GO TO 24 

23 PA=PC( I ) - ( PC ( I )-PB)*EXP ( -ET A ) 

AMA= SORT (5 .*( ( l.+AMB**2/5 . ) * ( P B/P A ) ** . 2 85 7 1 - 1 . ) ) 

XANS(I)=(X( I) +X( 1-1 ) 1/2 . 

E TAN S=E T A/2 . 

PANS (I )=PC( I )-(PC( I )-PB)*EXP(-ETANS) 

AMAN S( I ) =SQKT(5.*( ( 1.* A MB **2/ 5 . )*(PB/P ANS ( I) )**. 28571-1 . I ) 
DPDSA=(PC( I )-PA)/( PC( I )-PB) *DP DS 

24 CP( I )=(PANS( I >-l.)/(AM**2*.7) 

C DWS ( I ) =CDW S ( 1-1 )+CP( I )*( YSAVE ( I ) **2-YSAVE< I - 1 ) ** 2 ) *3 . 1 41 6 
AMU= ATAN ( (1 . /A MANS ( I ) I/SORT ( 1 .-(1 ./AM ANS ( I ) )**2 ) ) +OEL ( I ) 

C=SQRT( (6.* PANS l I ) +1 . ) / < 7 . *AM**2 ) ) 

THE T A= AT AN < C/SORT ( 1 ,-C**2 I ) 

YS( I ) = IY(I )+(XS( I — 1 ) — X ( I) )*SIN(AMU)/COS(AMU)-SIN( AMU ) /CO S ( AMU ) *COS 
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1 (THETA) /SIN ( THETA ) *YS( I -1 1 ) / ( 1 S INC AMU >/COS( AMU) *COS ( THETA ) IS IN (T 
2HETA ) ) 

ANGLE < I >=CDELC I) + AS INC 1 ./AMANS < I ) > ) *57. 29577 
XSCI) = XSC I-l I+CYSC II-YSC 1-1 ) )/S IN (THETA >*COSC THETA) 

25 WRITE <6,40) XANS ( I) , P ANS < I ) , AMANS ( I ) , CP < I ) , XS < I ) , YS < I ) , ANGL E ( I ) , C 
1 DWS ( I ) 

CALL BOUND 
ICODE=ICOp£+l 

IF ( ICODE.GT .ITER) GO TO 27 
I =N 

DO 26 KK=1 r I 
XNE W = X ( KK) /12. 

CALL FTLUP C XNEW , DLSTNW , +1 , I . XT, DLST AR ) 

0LSTNW=DLSTNW*12. 

Y( KK )=YSAVE(KK)*SQRT( 1 . +2 . *DL S TNW/ YS AV E ( KK ) ) 

XOL ( KK) =XCKK)/XCN) 

26 YCJL <KK)=Y(KK)/XtN) 

GO TO 5 

27 DO 28 1=1, N 

28 Y( I ) =Y S A Vt ( I ) 

GO TO 1 
C 

29 FORMAT (12A6) 

30 FORMAT (lbH ERROR IN F AL G/ 2 X , 6E 1 7 . 8 ) 

31 FORMAT (17H IMAGINARY ROOT S , 3E 1 7 . 8 / 6 E 1 7 . 8 ) 

32 FORMAT (17H TWU EQUAL ROOT S/ 2 X , 6 E 1 7 . 8 ) 

33 FORMAT ( 1H 1 , // / 1 X , 12 A6 ) 

34 FORMAT (1X25H1ST ORDER SHOCK EXPANSION) 

35 FORMAT (III 32X 14H INPUT D AT A , 1 CX , 12 H IT ER AT ION I 2/// 37X2HM=F 5 . 

12/36X3HRN=E13.8//) 

36 FORMAT ( 1 2 XI 9H ST AR T I NG CONE - M=F 8 . 5 , 4X6 HTHE T A= F 8 . 3 , 4X7HDEL SOR= F 8 

i .Min 

37 FORMAT (24X1HX,14X1HY» 1 3X5H DELTA//) 

38 FORMAT ( 2\. XF 8 . 4 , 7X F 8 . 4 , 7 XF 8 . 4 ) 

39 FORMAT ( / / / / 36 X6 HJ UTPUT / / 1 5 X 1 7 HI NV I SC I D SOLUT ION/ / 1 1 X4HX , I N , 5X8HP/ 

1 P( INF ) , 5X1HM»1GX2HCP»6X7H X- SHOCK, 5X7HY-SHOCK,1X20HCHARACTER I ST IC A 
2NGLE *2X9 HO /Q, SQ IN) 

40 FORMAT <7X,F8.4,3X,F9.5,3X,F7.4,3X,F7.4,2X,F10.4,2X,F10.4,8X,F8.4, 

16X.F10.6) 

END 


SUBROUTINE BOUND 

DIMENSION THE TAT 1200), DTHODX ( 200 ) , DMODX ( 200 ) , CFTC 200), SUM C 4 ) , 
1F0FXC4), 1HET2C20Q), CFT2I200), DTHDX2 C 200 ) » XTC200), DYODXC2CO), 
2TH0DLT (20o) , XC200 ), YSAVEC230 ) , PANS ( 203 ) , XANSC20D), AM AN S ( 200 ) , 

3 MT2C200), DLSTARC 200), YC200), VC 200), DELTT ( 200 ) , DEL C 200 ) , CDWS 

4 C 200 ) , XTTC2G0), ANGLE ( 2D0 ) , XLC200) 

COMMON XANS, AM AN S , MT2 , DL S T AR , N , XT , S , CF T , CFT 2 , V , DE LS OR , AM,RN, THODLT 
1 ,X, I tDELTI , Y, 0,YSAVE, PANS, DEL,T, CDWS, XL, THET AT , XT T , ANGL E ,CP C 200 ) , X 
20L ( 2 00 ) , YOL (200) ,Y MAX, ITER, [CODE , ST ( 200 ) 

REAL MT , MI 2 , NUTUAT , MP1 , MM1 , M , NUT AT 2 
EXTERNAL FUNC 
XT ( 1 )=XANSC 1 )/ 12. 

CFT2 Cl) =0. 

CFT ( 1)=0. 

XTT ( 1) =0. 

VC 1 ) =YS AVE ( 1 )/24. 
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00 1 1=2, N 

XT { I I =XANSI I )/12. 

CFT2 ( I)=0. 

CFTI I ) =0 . 

XTT( I J = XT( I )+{ CELT T ( I) -DLSTAR( I ) ) *SQRT I AM ANS ( I )**2-l. ) 

IF ( XTT( I } .GT.XTTI I- 1 ) ) GO TO 1 
XTT( I ) =XTT( 1-1 )+X(N)/12C‘O0. 

VI I )=< YSAVEl I) +YSAVEI l-l ) )/24. 

NUT AT2=AM*( i. + AM** 2/5. ) ** < -2 . 5 I /RN* ( 1.+216./T )/Il.+216./T/( l.+AM** 
12/5. I ) 

NUT0AT=NUTAT2*EXP< OELSOR) 

DO 2 1=1 ,N 

XL< I >=XT< I )/X(N)*12. 

A1=0 . 

CALL FTLUP ( X ANS ( I » , SS , + 1 , N , X , ST ) 

B1 = SS/12 . 

CALL MG A USS I A1 , B1 , 1 5 , SUM , F UNC ,FOF X , 4 ) 

THE T 2 ( I)=(l.+.4/2.*AM**2)**<-.2>*{ . 0076*N UT AT 2**1 ,167)*(SS/12.))** 
1 ( .857) /( AM** ( . 145) ) 

CALL FTLUP ( XT ( I ) , AMANX , + 1 , N , XTT , AM ANS ) 

THE TATI I ) = ( 1 . + .4/2 .*AMANX**2) **3/ I AM AN X** 3* VI I ) >*( . 00 76*NUT0 AT** . 1 
167/1 I 1. +.4/2.*AMANX**2 )**(-. 262>>*SUM( 1 ) )**. 857 
00 3 1 = 1 ,N 
XP 1 = XT ( I ) + . L C 1 
XM 1 = X T ( I ) - . Oul 

CALL FTLUP I XP 1 , TH ETP 1 , +1 , N , XT , THET AT ) 

CALL FTLUP ( XM 1 , TH E T M 1 , + 1 , N , X T , THE T A T ) 

CALL FTLUP ( X P 1 , MP 1 , +1 , N , XT T , AM AN S ) 

CALL FTLUP { XM1 , MM 1 , +1 , N , XT T , AMANS ) 

CALL FTLUP I XP 1 , TH ET P2 , + 1 , N , X T , THE T2 ) 

CALL FTLUP (XMi ,TH£T M2, +1»N»XT ,THET2 ) 

CALL FTLUP I XP1 , YP 1 , +1 , N , XT , V ) 

CALL FTLUP I XMI , YM1 , +1 , N , XT , V ) 

XD=XT< 11*12. 

CALL FTLUP ( XD,DL, +l,N,X,OEL) 

DYOOXI I)=(YP1-YM1) /(XP1-XM1 )*COS(OL) 

OTHOOXI I )=( THETP1-THETM1 ) / { XP 1-XM1 )*COS( DL ) 

0THDX2I I > = < THET P2 -THET M2 )/( XP1-XM1 )*COS(OL) 

OMOOXI I ) = ( MP 1 -MM 1 ) /( XP1-XM1 )*COS( DL ) 

OMOOXI 1 ) =E)MGDX I 3 ) 

OMOOXI 2 ) =OMUOX( 3) 

DO 1C 1=1 ,N 

CALL FTLUP ( XT II) , M , + 1 , N , X T T, A MANS ) 

CALL FTLUP I XT ( I ) , THET T , +1 , N , X T , T HE T AT ) 

A= 1 ./( l.+.2*M**2) 

6=1 .-A 

NSUM 1 = ( 0-1 . J /2 . 

SUM 1 . 

DU 4 I I = 1, NSUM 1 

SUMl=SUMl + ( B**I I/I 2 . *FLO AT I I I ) **2 + F L 0 A T I t I ) ) ) 

NSUM2=0- 1 . 

SUM 2 = 1. 

DO 5 11=1, NSUM2 

SUM2=SUM2+ I B**I I/I I FLOAT! II ) + 1 . ) * I 2 . *FL OA T I II ) + 1 . ) ) ) 

CALL FTLUP I M , THUD L , +1 , 5 1 , MT2 , THUOL T ) 

DEL T T ( I ) =1./TH0DL*THETAT( I ) 

DEL0R=DELTT I I ) /VII) 

CTH0DL=0/2 . *A/B**( I 0+ 1 . ) /2 . ) * I 2 . -AL OGI A ) + 1 . / S QR T! B ) * AL OG I I 1 . - SOR TI 
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1 B ) ) / ( 1 . 4 - SORT (B) 1 l-SUMl I +DEL0R*0/2.*A/B**0*( 1 ,/B*ALOG< A )-l./SQRT( B) 
2*ALOG( ( 1 .-SQRT(B) ) / ( 1 . + SORT ( B ) ) J-SUM2) 

IF ( (1.-(THOOL/CTHODl)> .LT. . 00 1 » GO TO 7 

THODL=CTHODL 

GO TO 6 

7 SUM3=0. 

DO 8 I I = 1 » NSUM1 

8 SUM 3 = SUM 3+ B** I l/FLOATI l I ) 

SUM4=1 . 

DO 9 11=1 *NSUM2 

9 SUM4=SUM4+8<‘*I 1/(2 ,*FLOAT( I I) +1. ) 

DSTODL = l . + 1 . /<l.*DEL0K4-0/2.*4/ ( B** ( (0 + 1,)/ 2.) ) *{ SUM3 + AL OG ( A ) ) +DELOR 
1*0*A/B**0*( SUM4-1. / (2. *SQRT ( B) ) *ALOG( ( 1 ,+SQRT <B) ) / ( 1 .-SQRT ( B ) ) ) ) 
DFLSTR=DSTODL/CTHODL 
DLSTAR( I )=DELSTR*TH£TT 
CFT2 ( I >=2.*DTHDX2( I ) 

10 CFT ( I I =( (DTHOUKt l » +THETT*( 1 ,/V ( I I *DYODX( l I +( 2 ,-M**2«-DELSTR) / { M* ( 1 . 
l+.4/2.*M**2 ) )*OMODX( I) ) 1*2. )*( M/AM)#=F2*PANS( I > 

A1 =0 . 

Bl = X ( N ) / 12 • 

CALL MGAUSS ( A 1 , B1 , 1 5 , S UM , FUNC , FOF X , 4 ) 

CFF I N=SUM ( 2 ) /SUM( A ) 

CFF I N2 = S UM ( 3 )/SUM(4J 
CDW=CDWS(N)/S/144. 

CDFR =SUM ( 2 »/S 
CDTOT=CDW+COFR 

WRITE (6 t il ) (XL( I ) ,THETAT( I) , CFT ( It ,CFT2< I) ,DLSTAR( I) , DELTT( I ) , 1 = 
11, N) 

WRITE (6,12) CFFIN.CFFIN2 
WRITE (6,13) COW, CDFR, CDTOT 
RETURN 
C 

11 FORMAT ( /21X16HV ISCUUS SOL UT I 0 N/ / 1 6X 3H X /L , 7 X 8 HTHE TA , FT3 X46HB00 Y-CF 
1, LOCAL FLAT PLATE CF, LOCAL DEL ST AR , F TAX 8 HD EL T A , FT/ ( F 2 1 . 5 , F 1 3 . 8 , F 
214.8,2F17.8,F14.8) ) 

12 FORMAT ( /2AXLAHBUDY-CF, TOT A L= F 1 1 . 8 , 7 X20HFL AT- PL A T E C F , TO TAL = F 1 1 . 8 I 

13 FORMAT ( /lBXBHCD,wAVt=Fll.8,7X12HCD,FRICTlON=Fll . 8 , 7 X 9HC D, TO T AL = F 1 
1 I .8 ) 

FND 


SUBROUTINE FUNC (XX,FOFX) 

DIMENSION XANS ( 203 ) , AMANS(20C), MT2(2C0), DLSTAR(290), XT(20C), C 
IF T( 2 00 ) , CFT 2 ( 20u ) , V( 20G > , PANS 1209 ), F0FXI4), Y(200), YSAVE ( 200 ) 
2, THODLT (2uC) , X(20u>, DELTT ( 200 ) , DEL(200), CDWS(230), XTT(2G0), 

3 ANGL E ( 2 00 ) , XL(2ui>), THET AT ( 200 ) 

COMMON X ANS , AMANS , MT 2 , DLST AR,N , XT , S »CFT ,CFT2 , V , DELS OR, AM,RN, THODLT 

1,X,I , DELTT, Y,0, YSAVE, PANS, DEL, T, COWS, XL. THETAT, XT T, ANGLE, CP (2C0),X 

20L ( 2 00 ) , YUL l 2c0 ) , Y MAX, [TER, ICODE , ST ( 200 ) 

REAL MM 

CALL FTLUP ( XX , Y VA L , + 1 , N , X T , V I 
IF (YVAL.LT.G) YVAL=G. 

CALL FTLUP ( XX , MM , +1 , N , X T , AM ANS I 
CALL FTLUP ( XX , CF , + 1 , N , X T , C FT ) 

XD=XX*12 . 
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CALL FTLUP ( XD ,DL, +1 * N , X , DEL ) 

FDF X ( 1 ) =(MM**3.33* YVAL**1 ,167/ ( { 1 . + . 4/ 2 . *MM**2 ) **4 ) J /COS (OL ) 

FOFX (2»=2.*3 .1416*YVAL*CF 

CALL FTLUP ( XX ,CF2 , +1 , N , XT, CFT 2 ) 

FOFX (3)=2.*3.1416*YVAL*CF2 
FOFX (4) =2.*3.1416*YVAL 
RETURN 
END 


SUBROUTINE BILUP ( TAB I , T AB J , T A B I J , T AC I J , T AD I J , T AE I J , N l , NJ , V AL I , V AL 
lJtBVALl *CVAL1»DVAL1»EVAL1 I 

C A TWO DIMENSIONAL TABLE LOOK-UP FOR TWO VARIABLES. 

C INPUT TABLES ARE T AB I J ( I , J ) , T AC I J ( l , J ) , T AD I J ( I , J ) AND TAEIJ(I,J) 

C AS FUNCTIUNS OF TABI(I) AND TABJ(J). THE TWO DEPENDENT VARIABLES A 

C LINEARLY INTERPOLATED SIMULTANEOUSLY FOR INPUT VALUES OF VALI AND 

C V AL J RESULTING IN ANSWERS BVAL1, CVAL1, OVAL1, EVAH. ERROR 

C SIGNALS ARE GENERATED WHEN THE TABJ(J) TABLE IS EXTRAPOLATED. 

DIMENSION TABJ(ll) , TABI(13,11>, TAB I J ( 13 » II ) , TAC I J ( I 3 , II ) , TADIJ 
1(13,11), TAE IJ( 13, 11) , TB I SL ( 2 ) » TBIJK2), TCIJK2), TDIJK2), T fc I 
2J1 ( 2 ) 

TB J S L=0 . 0 
KK=2 

IF ( TABJ( 1 ) .LT .TAB J(2) ) GO TO 4 
DO 2 J= 1 ,N J 

IF ( VAL J-T AB J ( J ) ) 2,1,3 

1 TBJSL=1. .0 
GO TO 8 

2 CONTINUE 

IF (J.EQ.NJ) WRITE (6,22) 

3 IF ( J.GT .1 ) GU TO 8 

WRITE (6,23) 

J=2 

GO TO 8 

4 DO 6 J=1 »NJ 

IF ( VAL J-T AB J ( J ) ) 7,3,6 

5 TBJSL=1 .0 
GO TO a 

6 CONTINUE 

IF (J.EQ.NJ) WRITE (6,22) 

7 IF (J.GT. 1) GU TO 8 

WRITE (6,23) 

J = 2 

8 IF ( TABI (1 , J) . GT.T ABI ( 2, J) ) GO TO 12 
IF ( VALI .LT.TABKl , J) ) GO TO 10 

IF ( VALI .GT.TABI ( N I , J) ) GO TO 11 
DO 9 1=1, INI 

IF ( VAL I — T A B I ( I,J) ) 17,16,9 

9 CONTINUE 

10 1=2 

GO TO 17 

11 I =N I 

GO TO 17 

12 IF ( VALI .GT.TABK 1 , J) ) GO TO 14 
IF ( VALI.LT.TABKNI, Jl I GO TO 15 
DO 13 1=1, NI 

IF ( VALI-TABI ( I , J ) ) 13,16,17 

13 CONTINUE 
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14 1=2 

GO TO 17 

15 I=N I 

GO TO 17 

16 TBI SL( KK ) =0 • 0 

TBIJ1(KK)=TAB1J(I , J) 

TCIJ1(KK)=TACI J(I, J) 

TD I J 1 { KK ) =T AOl J( I, J) 

TEI J1IKK )=TAE1 J( I , Jl 
GO TO 18 

17 TBI SL(KK» = { VALI-TABK 1-1, J) >/(TABl(I , J I -T ABI ( 1-1 , J I) 

TBI J1(KK)=TBISL(KK)*(T ABIJI I*J )— T AB IJ{ 1-1 , J M +T AB 1 J ( I- 1 , J ) 
TCI JKKK )=TBISL(KK)*(TACI J< 1 , J I-TACI J( 1-1, J) ) +TACIJI I- 1 * J ) 
TDIJ1(KK) = TBISL(KK)*(TADIJ( I,JI-TADIJC I -1 , J ) > +T AD I J { I- l, J) 
TEIJ1(KK)=TB 1 SL ( KK I *{ T AE I J ( I , J )-T AE I J ( I -1 » J > ) +T AE I J < I -1 ♦ J ) 

18 IF (T8JSL.EQ.G.0) GO TO 19 
B VAL 1=TB I Jl ( KK I 
CVAL1=TCIJ1 ( KK ) 

DVAL1=TDIJ1( KK I 
EVAL1=TEIJ1(KK) 

GO TO 21 

19 IF (KK.EQ.l) GO TO 20 
KK=K K- 1 

J = J-1 
GO TO 8 

20 J = J+1 

TBJSL=( VALJ— TABJ(J-1 ) ) /(TABJ (J )-TABJ<J-l) ) 

BVAL1=TBJSL*(TBI Jl I2I-TBI Jl ( 1 ) I tTBI Jl (11 
CVAL1 = TBJSL*UCI Jl <2>-TCIJl(in +TC I J 1 ( 1 ) 

DVAL1=TBJSL*(TDIJ1 (2J-TOIJ1 ( 1 I H-TDIJK 1 ) 

EVAL1=TBJSL*(TEI Jl (2J-TEI Jl( 1 ) I+TEIJK 1 I 

21 RETURN 

C 

22 FORMAT I/20X20HHIGH J EXTRAPOLATION) 

23 FORMAT I/20X19HLUW J EXTRAPOLATION) 

END 


SUBROUTINE MGAUSS ( A , B , N , SUM , F UNC , F OFX . NUMUtR ) 

C THIS SUBROUTINE INTEGRATES FROM ZERO TO ONE 

DIMENSION U ( 5 ) . R(5). SUM ( 1 ) , FOFX(l) 

DO 1 LL= 1 , NUMBER 
1 SUM C LL ) = 0 • 0 

IF (A.EQ.B) RETURN 
U(1 ) = . 425562830509184 
U ( 2 ) = .283302302935376 
U ( 3 ) = . 1 60295215850488 
U ( 4 ) = . 06746831 6656508 
U ( 5 ) = . 0 1 304673574 1414 
R ( 1 ) = . 1 47 7621 12 35 73 76 
R (2 )=. 1 34633359654998 
R < 3 ) = . 1 095431 8 1 25799 1 
R ( 4 ) = .074725674575290 
R ( 5 ) = .033335672 1 54344 
F I NE=N 

DELT A = F I NE/ (B-A ) 
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DO 3 K=1 ,N 
XI =K-1 

F I NE=A+X I /DELTA 
DO 2 11=1.5 
UU = U( I I ) /DELTA+F I NE 
CALL FUNC (UU.FOFX) 

DO 2 JOYBOY= 1 .NUMBER 

SUM ( JOYBOY ) =R C I I ) *F OFX (JOYdOY) + SUM ( JOY^OY) 
DO 3 JJ=1., 5 

UU= ( 1 . O-U ( JJ ) ) /DELTA+F I NE 
CALL FUNC (UU.FOFX) 

DO 3 NN=1 .NUMBER 

SUM(NN)=R(JJ)*FOFX(NN ) + SUM (NN ) 

DO 4 I JK=1» NUMBER 

SUM (IJK) = SUM (IJK) /DELTA 

RETURN 

END 


SUBROUTINE FALG ( C OEFF S . N , I . ROOT , TEMP , I ERR ) 
DIMENSION COEFFS(l). TEMP ( 1 ) , ROOT ( 1 ) 
DIMENSION XID1FFI2). RD IFF ( 2 ) , APPROX<3) 
COMPLEX F ,FPR, APPROX, TEMP, ROOT 
COMPLEX tempm 
COMPLEX RELTST 
NS AVE = N 
I ERR=0 
I B = 0 

I CLEAN = 2*N+2 

clear working area 

DO 1 LLL= 1 , I CLEAN 
TEMP (LLL ) =0. 0 

clear root storage 

DO 2 LLL = 1 , N 
ROOT (LLL ) = 0 • 0 

CONSTANTS to test 
CONVERGENCE 
CONST = . 1 E-6 

OVERFLOW 
0 V CON= 1 .El 50 

MAGNITUDE OF ROOTS 
RCONST = 1 . E-2 1 

JONJON=0, FIRST ITERATION 
JON JON = 0 

CHECK CONSTANT TERM FOR ZERO 

JJ J= 1 
NCO=N+ 1 

IF ( I .NE. 1 ) GO TO 4 

COMPLEX COEFFICIENTS 
NCO=2*NCO 

IF ( COEFFS ( NCO- 1 ) • NE . 0 . ) GO TO 5 
HERE IF REAL COEFFICIENTS 
IF (COEFFS (NCO ) .NE.O. ) GO TO 5 
ROOT = ZERO 
ROOT ( JJ J ) =0. 

NCO=NCO- 1 
JJ J=J JJ+ 1 
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C REDUCE DEGREE AND IF 1 , STORE ROOT AND EXIT 

N = N— 1 

IF (N.NE. 1 ) GO TO 3 
ROOT ( JJJ ) =0. 

GO TO 38 
C 

ENTRY FIRST AND SECOND ITERATIONS 

J = JJ J 

NTERMS=N+ 1 
KCON J = 0 

clear approx 
approx ( 1 > = o . o 

APPROX ( 2 ) =0 . 0 
APPROX ( 3 ) =0 . 0 
IF (I.EQ.l ) GO TO 7 

REAL COEFFICIENTS 
DO 6 I JFF=1 ,NTERMS 

TEMP ( I JFF ) = CMPL X ( COEFFS ( I JFF ) , 0 .0 ) 

GO TO 9 

COMPLEX COEFFICIENTS 
DO 8 I I I X= 1 i NTERMS 

TEMP ( I I I X ) =CMPLX < COEFFS ( 2* I I I X- 1 ) , COEFFS (2* 1 I I X ) ) 

CHECK LEADING COEFFICIENT FOR 0 OR 1 
TEMPL = REAL ( TEMP ( 1 ) ) 

IF (TEMPL.NS.O. ) GO TO 10 

IF REAL IS ZERO, CHECK IMAGINARY 
TEMPL = A I MAG ( TEMP ( 1 ) ) 

IF CTEMPL.NE.O* ) GO TO 10 

LEADING COEFFICIENT ZERO IERR=1 

I ERR= 1 
GO TO 38 

DIVIDE BY LEADING COEFFICIENT 

0 TEMPM=TEMP ( 1 ) 

DO 1 1 LL A = I , NTERMS 

1 TEMP (LLA ) =TEMP (LLA J/TEMPM 

KCON J= 1 , TRIAL VALUE=CON JUGATE 

2 IF (KCONJ.NE.O) GO TO 13 

FIRST TRIAL VALUE 
APPROX (1 ) = ( . 0 1 , . 0 1 ) 

D I FFERENT I ate 

3 DO 1 4 I I = 1 , NTERMS 
XPON = NTERMS— I I 
NNOW = I I +NTERMS 

4 TEMP (NNOW ) =XPON*TEMP ( I I ) 

NPON = NTERMS — 1 

KA=0 FOR FIRST TRIAL VALUE 

K A = 0 

JONJON= 1 SECOND ITERATION 

5 IF (JONJON.EQ. 1 ) APPROX t 1 )=ROOT(J) 

CLEAR RDIFF,XIDIFF 
CO 16 LLL = 1 , 2 
RD I FF (LLL ) =0 .0 
5 X I DI FF ( LLL ) = 0. 0 

ROOT EVALUATION 
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C MAXIMUM ITERATIONS =120 

17 L = 2 

PARTR1 =REAL ( APPROX ( 1 ) ) 

PARTM 1 = A I MAG (APPROX ( 1 ) ) 

DO 25 K=2, 121 
C EVALUATE F(X) 

18 F= (O.O.C.O > 

DO 19 11 = 1 , NTERMS 

F = APPROX (L-l ) #F + TEMP ( I I ) 

XF=ABS (REAL (F ) ) 

YF = A3S ( A I MAG (F ) ) 

C CHECK FOR OVERFLOW 

IF ( XF . GT . OVCON. OR. YF.GT.OVCON ) GO TO 26 

19 CONTINUE 

C EVALUATE FPRIME(X) 

FPR= (0.0,0 *0 ) 

CO 20 JJ= 1 , NPON 
NN0W = JJ+NTERMS 

FPR= APPROX (L-l ) *FPR+ TEMP (NNOW ) 

YFP = ABS ( A I MAG ( ^PR ) ) 

XFP = ABS (REAL (FPR) ) 

C CHECK FOR OVERFLOW 

IF ( XFP .GT .OVCON. OR. YFP.GT . OVCON ) GO TO 26 

20 CONTINUE 

C SEE IF FPR]ME=0 

IF ( XFP .EQ . 0 • 0 • AND. YFP .EQ. 0 • 0 ) GO TO 26 
C IF NOT ZERO, NEW APPROXIMATION 

APPROX ( L ) =APPROX (L- 1 l-F/FPR 
PARTR2 = REAL ( APPROX ( L ) ) 

PARTM2 = A I MAG ( APPROX ( L ) ) 

C SET EITHER PART TO ZERO IF LESS THAN l.E-21 

IF ( ABS (PAPTR2 ) .LE.RCONST ) PARTR2=0. 

IF (ABS (PARTM2 ) .LE.RCONST ) PARTM2=0. 

IF ( PARTR2 • EQ • 0 . . AND .PARTM2.EO.O. ) GO TO 21 
GO TO 22 

ZERO ROOT 

1 IF (L.EQ.3) APPROX ( 2 ) =APPROX (3 ) 

GO TO 30 
C 

22 RD I FF ( L- 1 ) =ABS ( PARTR2-PARTR 1 ) 

XIDIFF(L-1 )=ABS (PARTM2— PARTM1 ) 

IF (L.EQ.3) GO TO 23 

L = 3 

PARTR1 =PARTR2 
PARTM1 = P ARTM2 
GO TO 18 
C 

C TEST 1 

23 IF ( (RD IFF (2 >+X ID IFF <2 ) ) .LT. (RD IFF < 1 )+X ID IFF < 1 ) ) ) GO TO 24 

C TEST 2 

RELTST = (APPROX ( 3 ) -APPROX (2 ) ) /APPROX ( 3 ) 

DIFFR = ABS (REAL (RELTST ) ) 

D I FFX I = ABS ( A I MAG (RELTST ) ) 

IF (DIFFR. LT, CONST. AND. DIFFXI .LT. CONST ) GO TO 30 

24 APPROX ( 2 ) =CMPLX (PARTR2,PARTM2 ) 

PARTR 1 =PARTR2 

PARTM I =PARTM2 
RDIFFf 1 ) = RD I FF (? ) 

25 X IDIFF ( 1 ) =XI OIFF (2 ) 
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c 

c 

c 

c 

c 

26 


C 

27 
C 

c 

c 

c 

28 
C 


c 

29 

C 

C 

30 
C 

C 


31 

C 

C 

c 


32 


C 


C 


C 

33 

C 


MAXIMUM ITERATIONS EXCEEDED OR 

OVERFLOW OR 

FPRIME=0 

TRY AGAIN WITH SECOND TRIAL VALUE 
IF ( JONJON.EQ. 1 ) GO TO 27 
IF (KA.E0.i05) GO TO 28 
APPROX ( 1 ) = ( 1 . , 1 . ) 

KA=1 05 
GO TO 17 

SECOND ITERATION NONCONVERGENT ROOT IERR=3 
I ERR=3 

STORE RESULT AND IMPROVE NEXT ROOT 
GO TO 34 

FIRST ITERATION ROOT R NONCONVERGENT IERR=2 
IMPROVE ( R— 1 ) ROOTS 
I ERR=2 

I B =LAST CONVERGENT ROOT 

I B = J 

ROOT ( J ) = APPROX (2 ) 

IF ( IB.NE. 1 ) GO TO 29 

FIRST ROOT FAILED RETURN 

GO TO 38 
JON JON = 1 
GO TO 5 

STORE ROOTS 

IF (JONJON.EQ. 1 ) GO TO 34 

HERE IF FIRST ITERATION 
ROOT ( J )= APPROX ( 2 ) 

REDUCE POLYNOMIAL BY SYNTHETIC DIVISION 
NTERMS = NTERMS— 1 
DO 31 I K=2 i NTERMS 

TEMP ( I K ) = ROOT ( J ) * TEMP I I K- 1 ) +TEMP (IK) 

CONT I NUE 

NEXT ROOT IF COMPLEX COEFFICIENTS 
IF ( I . EQ. 1 > GO TO 33 

HERE IF REAL COEFFICIENTS 
IF (KCONJ.LQ.O) GO TO 32 

RESET KCONJ IF ROOT IS CONJUGATE OF PREVIOUS ROOT 
KCON J = 0 
GO TO. 33 
x = real ( root < j > > 

Y = A I MAG ( ROOT ( J ) ) 

IF ( X • EQ » 0 • ) GO TO 33 

SEE IF REAL OR COMPLEX 
IF (ABS ( Y/X ) .LE . 1 .E-l 0 ) GO TO 33 

COMPLEX ROOT TRIAL V ALUE= CON JUGA TE 

APPROX ( 1 ) =CONJG (ROOT ( J ) ) 

KCON J= 1 

NEXT ROOT 

J = J+ 1 

IF (J.NE.NSAVE) GO TO 12 
LAST ROOT 

ROOT ( J ) = — TEMP ( 2 ) / TEMP ( I ) 

JON JON = 1 
GO TO 5 
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c 

C IMPROVED ROOT FROM SECOND ITERATION 

34 ROOT ( J ) =APPROX ( 2 ) 

X=ABS (REAL (ROOT ( J ) ) ) 

Y = ABS ( A I MAG (ROOT ( J ) ) ) 

C SET REAL OR I MAG. TO ZERO IF LESS THAN RCONST 

IF (X. LT. RCONST. AND. Y.LT. RCONST ) GO TO 35 
IF (X.LT .RCONST ) FOOT ( J )=CMPLX (0 .0 i A I MA6 (ROOT ( J ! ) ) 
IF (Y.LT. RCONST ) ROOT ( J ) =CMPLX ( REAL ( ROOT ( J )), 0 . 0 > 
GO TO 37 

35 IF (X.GE.Y) GO TO 36 

ROOT ( J ) =CMPLX ( 0 .0 . A 1 MAG (ROOT ( J ) ) ) 

GO TO 37 

36 ROOT ( J ) =CMPLX ( REAL ( ROOT ( J ) > , 0. 0 ) 

C 

C IF(R-l) ROOTS I MPROVED .RETURN 

37 IF (J.EQ.Id) GO TO 38 

J = U+I 

C IF N ROOTS I MPROVED . RETURN 

IF (J.LE.NSAVE) GO TO 15 

38 N = NS A VE 
RETURN 
END 


SUBROUTINE FILUP (X,Y ,M,N, VAR! ,VARD) 

01 MELS I UN VABI(i), VARO(I), V<3), YY(2) 

IF ( M . F 0 . m Ai'jO . N . E y . ^ J o ) TO 14 

IF ( M » F; 0 1 0 , 0 . I'm . Nc .0 ) GO TO 21 

IF ( N . L F • I A B 3 ( H ) ) GJ lU 21 
IF (M.GT.3) GU TO 8 
C M.LT.-j 

00 1 I Y Y = 1 , N 

1 = I YY 

IF (VAR T ( [ I— XI 2 , 2 o , 1 

1 CONTINUE 
[ = f\l f M 

C IF X.LT.x(N) fcXfKAPJLATt 

IF (M.FW.-l ) Gu I U 3 
GO TO 7 

C IF Y.GT.M 1) , EXI RAPOLATF 

2 IF ( I . EQ.l . A.Tbi.i-l.Lg.-iI GO TO 3 
IF (I .FQ.i . AimJ.M. tw.-2) GO TO 7 
IF IM.Nt: .-1 ) Gu TU 4 

C M = -l 

1 = 1-1 

3 IF (VAR I ( II .Lt.VAKK 1+i) ) GO TO 21 
GO TO 12 

M= -2 

IF (I.Nr.M) bu TU t) 

I = N-2 

Gil TO 7 

C C 0 M P A K F n' I I H NeXT 

5 I F ( VAR II I ■*- 1 ) — X ) o.Xl.zl 

6 1=1-1 

IF ( I , E 0 . 1 ) GJ TU 7 
C SEE WHICH THKlL 

IF ( <VARI( I -i I -X ) .LT. (X-VARI ( I +2 > )) 1=1-1 

7 IF ( VAR I (I) ,U . VARI ( 1+1 ) ,!JR .VAR I ( 1+1 > . LF. VAR I ( 1 + 2) ) GO TO 21 

GO TO 1 7 
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C M.GT.U 

8 00 9 I YY = 1 ,N 
I = I YY 

IF ( X— VAR l(I)| 10 j| 9 

9 CONTINUE 
I = N-M 

C IF X.GT.Xdv) »EXTKAP01ATE 

IF IM.EO.U GO TO a 
GO TO 16 

C IF X, LT. XI 1) .EXTRAPOLATE 

10 r F I I .60. 1 . AnD.M. tu.i ) GO TO 11 
IF ( I .EO.l .A nu.M.cU.c 1 GO TO 16 
IF < N.NC.U GO Tu Li 

C M= 1 

1 = 1-1 

11 IF (VAR I { 1+1) . L £ . V AR 1 (1) ) GO TO 21 

C LINEAR 

12 Y=(VARD{II#(VARH I+1)-X)-VARD( I +1 )*{ VAR I( D-K) I /( V4R !( PI I -VAR I ( I ) 
1) 

RETURN 

C M= 2 

13 IF (I.NE.N) Gu TU 14 
I = N-2 

GU TO lo 

C COMPARE With NEXT 

14 IF (X-VARH I fin 15,*:i,*:l 

15 1=1-1 

IF (I.cO.ll bU TO lb 
C SET WHICH TEiRct 

IF UX-VARI (I-i )l,Ll. I VARl (!+2)-XM 1 = 1-1 

16 IF ( VAR I ( I +1 } .LE.VARI (1 ) .OR. VARI ( 1+2 I ,LE. VAR I ( 1+1 )) GO TO 21 

C S ECONO UKDlR 

17 VI 1 ) = V AR I ( I )-X 

V( 2 ) = V A R I ( I+i ) — X 
VI 3 ) =V A R I ( 1+2 ) - X 
K= I 

00 IP. J = l,2 

YY I J ( = ( VARiM M* V( J+X) -VAR0(K+1 >*V < J ) )/< VAR I (R + l ) -V AR I (K > ) 

18 K = K + 1 

Y= ( YY( 1 )#V< 3 )-Y Y( *; ) *V ( 1 )) / ( VAR I ( 1 + 2I-VAP I ( l I } 

RETURN 

C ZERO URUtk <Y= Y( 1 ) ) 

19 Y=VAPD(i) 

RETURN 

C Y = Y(U 

20 Y = VARD( I I 
RETURN 

C 

C ERROR PRINT 

21 PRINT 22 
PRINT 23, M,N,X 
IF (N.EQ.O) si UP 
IF IM.EQ.O) STUP 
PRINT 24 

STOP 

C 

22 F0KM4T (/31H eRKuK WAS ENCOUNTERED IN FTLUP) 

23 FORMAT { l X , 2iH = 1 6 ,6 X , 2HN= I 5 , 5X , 2HX =E 20 . 8 ) 

24 FURMAT (19H I ABL c. UUT UF ORDER) 

END 
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Description of Input Data 

A single case consists of the determination of the inviscid flow field, body surface 
conditions, and boundary-layer growth along a meridian line from the nose to the base of 
the body of revolution. The Mach number, unit Reynolds number, and free-stream stag- 
nation temperature must be put into the program in order to define the free-stream flow 
conditions. It is necessary to use as input, in addition to the body geometry, a table of 
cone surface pressures, cone surface Mach numbers, cone shock angles, and cone entropy 
changes for a range of Mach numbers and cone half-angles. In addition, a table of the 
two-dimensional values of 8/6 for a range of Mach numbers is required for the desired 
velocity profile index. The number of iterations to be made between the viscid and 
inviscid solutions must also be specified. For the loading routine used in the program, 
any column except the first may be used on the input cards unless otherwise specified, 
and a decimal format is used for the input quantities unless a fixed point number is spec- 
ified. The FORTRAN names and a description of the required inputs are as follows: 


FORTRAN 

Available 
card columns 

Description 

WHAT 

1 to 72 

Identification card (any identifying information may be 
written on this card and will appear at the start of 
output for each case) 

$NUM 

2 to 5 

Arbitrary name required by the loading routine to define 
the block of input data 

AM 

2 to 80 

Free-stream Mach number 

RN 

2 to 80 

Unit Reynold's number, 12/model units 

N 

2 to 80 

Number of body coordinates specified (Fixed point number; 
200 points maximum) 

X(l) 

2 to 80 

Body x- coordinate array, model units 

Y(l) 

2 to 80 

Body y-coordinate array, model units 

TO 

2 to 80 

Stagnation temperature, °R 

O 

2 to 80 

Velocity profile index 

MT2(1) 

2 to 80 

Mach number array for two-dimensional boundary-layer 
parameter 9/8 (51 points) 

THODLT(l) 

2 to 80 

9/8 array (51 points) 
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FORTRAN 

Available 
card columns 

Description 

MTAB(l) 

2 to 80 

Mach number array for cone tables (11 points) 

DELTAT(1,1) 

2 to 80 

Cone half-angle array for cone tables (two-dimensional 
array; 13 points columnwise, 11 points rowwise), radians 

PCT(1,1) 

2 to 80 

Cone surface pressure Pc/p^ (two-dimensional array; 
13 points columnwise, 11 points rowwise) 

AMIT(1,1) 

2 to 80 

Cone surface Mach number (two-dimensional array; 
13 points columnwise, 11 points rowwise) 

THET1T(1,1) 

2 to 80 

Cone shock angle (two-dimensional array; 13 points 
columnwise, 11 points rowwise), radians 

DLSOR(l,l) 

2 to 80 

Entropy change across the cone shock, A?]/Rg (two- 
dimensional array; 13 points columnwise, 11 points 
rowwise) 

ITER 

2 to 80 

Number of iterations to be made between viscid and 
inviscid solutions (fixed point number) 

$ 

2 

Denotes the end of a case 


The system loading subroutine used in the program (NAMELIST) is quite flexible 
in that the order of the input cards is unimportant and successive cases can be run by 
repeating the WHAT and $NUM cards followed by only the changed parameters and a 
$ card. 
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Example Input Listing 

As an example of a set of input cards, the following is a listing of the inputs 
necessary to compute the viscid and inviscid flow parameters on the RM-10 body at 
M m = 2.2 and R M = 4.16 x 10® with two iterations specified: 


RM 10 BODY 
$NUM 

MTAB(l) -1.5,2. ,2. 5, 3. ,3. 5, 4. ,4. 5,5. ,6. ,8. ,10., 

DELTAT (1,11=0., .043 63 323,. 08726 645,. 13 08 99 67,. 17 45 3 29,. 21 81661,. 2 61 7993, 
.3054325, .3490658, .392699, .4363322, .479965, .5235987, 

0.,. 04363 323,. 087 26 6 45,. 13089 967,. 1745 3 29,. 21816 61,. 2617993, 

.3 0543 25,. 349065 8,. 3926 99,. 43 63322,. 479965,. 5 235 987, 

0.,. 0 43633 23,. 087266 45,. 1308 9967,. 1745 3 29,. 21 81661,. 2 617 993, 
.305 43 25,. 3490 658,. 39 26 99,. 43 63322,. 47 9965,. 5235987, 

0.,. 043 633 23,. 08726645,. 13089967,. 17 45 329,. 218 1661,. 2 617 993, 
.305 43 25,. 34 90 65 8,. 3926 99,. 4363322,. 47 9965,. 52 35 987, 

0.,. 0436 3323,. 08726645,. 13089 967,. 1745 329,. 21 81 6 61,. 2617993, 

.3 05 43 25,. 34 90 65 8,. 3926 99,. 4363 322,. 479965,. 5235 987, 

0 .,. 043 63323,. 087 26 645,. 13089 967,. 1745329,. 2 181 661,. 2617 99 3, 
.3054325, .3490 658,. 39 2699,. 43 63 322,. 47996 5,. 5235 987, 

0.,. 04 36 33 23,. 08726645,. 13 089967,. 1745 329,. 218 1661,. 2617993, 

.3 054325, .3 490 658, .3 92699, .43633 22, .479965, .5 235987, 

0., . 04363 3 23, .08 726645, . 13089967, .1745 329, . 2181661, .2617993, 
.305 43 25,. 349 0658,. 392699,. 4363 322,. 479 965,. 52 35 98 7, 

0 ., . 043 633 23, . 08726645, .13089967, .17 453 29, . 2181661, . 2617993, 
.305 43 25,. 3490 65 8,. 39 26 99,. 43 63322,. 479965,. 52 35 987, 

0., . 043 63 3 23, .087 26645, . 13089967, .17 45 329, . 2181661, .2617993, 
.3054325, .3490658, .392699, .4363322, .479965, .5235987. 

0 ., . 043 633 23, . 08726645, .13089967, . 17 45 3 29 , . 2 18 1661 , . 2 61 7 993 , 
.3054325,. 3 4906 5 8,. 39 26 99,. 43 63 322,. 4799 65,. 5235 987, 

PCT(1,1)=1,, 1.01 9 43 26, 1.06 24668, 1.1218836, 1.19498 23, 1.28 05 338, 1.378 04 93, 
1.48 75 334, 1. 60 94921, 1.7452173, 1.89 77323, 1 . ,1 ., 

1 ., 1.0 30 201 4, 1.095 0 48, 1.183 62 74, 1.29 2483 2, 1.4 20 2392, 1.5662 549, 
1.73013, 1.91 14873, 2. 1099574, 2. 325245, 2. 5572482, 2. 8063347, 

1.. 1.0 430054. 1.13381 27. 1.25 79479. 1.4118223. 1.59 43 90 8. 1.80 51 501, 
2.0435503,2.3088209,2.5999833,2.9158958,3.2553311,3.6170774, 

1.. 1.0575531. 1.17 7 91 35. 1.3434709. 1.5 51 093 6. 1.8 00 2508. 2. 090 53 85, 
2.4211572,2.7908532,3,1979721,3.6405103,4.1162105,4.6226939, 

1.. 1.07 368. 1.2270 21 1.1. 43 99168. 1.71 01 6 05. 2. 0377458. 2. 4222429, 
2.8624798,3.3565651,3.9019948,4.4957647,5.1344738,5.8144615, 

1.. 1.0 91 27 20. 1.2809 8 67. 1.547292. 1.88920 65. 2. 3 07 1757. 2. 8 00 6264, 

3. 3678 42 9, 4. 0061275, 4. 71 195 93, 5. 4811532, 6. 308 9994, 7. 1904022, 

1.. 1.110 2609. 1.33975 31. 1.665 6866. 2. 0884875. 2. 60891 72. 3. 2261054, 
3.9376627,4.7398834,5.6280810,6.5967119,7.6395491,8.7499419, 

1.. 1.13 058 48. 1.4 0331 5. 1.7 95 246 8. 2. 308 246. 2. 9 43 2 81 6. 3. 69 905 42, 
4.5722872,5.558163,6.6506071,7.8425504,9.1261264,10.492909, 

1.. 1.17 50571. 1.54 49088. 2. 0883718. 2. 81 0078 4. 3. 711 01 58. 4. 78 85 278, 
6.0371455,7.4492676,9.0156074,10.725597,12.5675,14.528819, 

1 ., 1.278 69 74, 1.8871 344, 2. 8137 638, 4. 0675451, 5. 6477433, 7. 5 47 1315, 

9. 7541022, 12. 253 764, 15. 028693, 18. 05 9352, 21. 3 2435, 24. 80 0093 6, 

1 ., 1.4015141, 2, 310 2 668, 3. 7291669, 5. 6 691 93 7, 8. 1254669, 11. 084206, 

14 .525 483, 18. 4 25 183, '22; 755 51 3, '27. 485-475, 37. 581175, 38. 00706, 
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AM1T(1 / 1)=1. 5, 1,4866884, 1.4579318,1. 4197316, 1.3748579, 1.3248992,1. 2707392, 

1. 2 127 162, 1.1506 488, 1. 0837 3 13, 1. 01007 03, l.,l., 

2. . 1.9808 61 5. 1.94 15 25 4. 1.8 91 2 335. 1.83 40453. 1.77 21 933. 1.70688 37, 

1.63 8 67 64, 1.56 77 572, 1.4 94 0884, 1.41 7468, 1.3375202, 1.2 53 5 953, 

2.5,2.4729472,2.4193976,2.3528628,2.2788768,2.2001554,2.1179422, 
2.0327606,1.9448503,1.8543543,1.7613748,1 .6659616,1 .5680718, 

3. . 2. 96276. 2. 891 38 32. 2. 8048038. 2. 71 01402. 2. 6103 929. 2. 5067 511, 
2.3997895,2.2899646,2.1777409,2.063583,1.9479054,1 .8310067, 

3.5,3.4500939,3.3571597,3.2467499,3.1275136,3.0025569,2.8731152, 
2.7400668, 2.6043443, 2.466913,2 .3286833, 2.1904152, 2.0526629 , 

4.. 3. 93477. 3. 81 645 43. 3. 6783916. 3. 5305 68. 3. 3761469. 3. 21 6 6972, 
3.0537628,2.8889908,2.7239783,2.5600783,2.3983186,2.2393739, 

4.5,4 .4166246,4 .269058,4.0994611,3 .9189007,3.730775,3.5374323, 
3.3414321,3.1453333,2.9513677,2.7612386,2.5760876,2.3965344, 

5. . 4. 895 5305. 4. 7 14 807 8. 4. 5 097079. 4. 29 21 847. 4. 0662668. 3. 835 593 2, 
3.6040616,3.3752078,3.1517935,2.935687,2.7279347,2.5289133, 

6. . 5. 84406 86. 5. 5852 47 4. 5. 2968787. 4. 99 2783 2. 4. 6801702. 4. 36 681 18, 

4. 0597 311, 3. 76 4078 5, 3. 4829712, 3. 217 8167, 2. 968793 4, 2. 735 25 8, 

8 .. 7. 7 01 8785. 7. 23 90804. 6. 7337 145. 6. 2103924. 5. 69226 72. 5. 198150, 

4. 73895 43, 4. 31 87 765, 3. 93 73 126, 3. 5918831,3. 2787 85 8, 2. 9 940481, 

10. ,9. 50464, 8. 77 25 147, 7. 9851813, 7. 19 77 58, 6. 456834, 5. 786 75 81, 5. 193 3727, 

4. 6723 71 7, 4. 2153623, 3. 81 31 01 9, 3 . 4569163, 3 . 13 91 94, 

THET1TC 1,1 > = . 7 297 92 25, .7 29 79225,. 73078309,. 73 489674,. 74 4660 72,. 76158205, 
.78592 36 5, .81 737375, .85571959, . 9013 2794, .9557 85 04, l.,l., 

.523 69 89,. 5236989,. 52524879,. 53140177,. 5446 48 27,. 5 65 18145,. 5919218, 

.6237 01 15,. 65966167,. 69 928055,. 7 42 32301,, 78882392,. 83913 45 3, 
.41168587,. 4 11 68587,. 414 25 491,. 42357796,. 44136238,. 46630615,. 4966 26, 

.5310 16 45,. 56864312,. 608 99456,. 65177 73,. 6968639,. 74427222, 

.34011 05 7,. 34011057,. 34409993,. 3 57 01 069,. 37 899374, .40760084, .44085123, 
.47752193,. 51687 237,. 5 5844619,. 601966 73,. 64728576,. 69436247, 
.2901713 4,. 2 90 17134,. 2 95 91943,. 31238763,. 33790128,. 3693 9419,. 4 0494077, 
.44344094,. 48 4237 53,. 5 26921 4 7,. 57123978,. 6170 4946,. 66429713, 
.25329402,. 2532 9 402,. 26104195,. 28078488, .30918 22,. 3 42 95689,. 3803207, 
.42028295,. 4 62 25228,. 5 05 85 84,. 550 87 169,. 5971618,. 64467684, 
.22495268,. 22 49 52 68,. 23483 71,. 25747371,. 2882 2452,. 3 2 383435,. 36266044, 

.40 38 0119,. 44671 75 4,. 4 91 07302,. 53666064,. 5833 64 5,. 63114, 
.2025178,. 202 51 78,. 21458047,. 2397 25 46,. 27 241 57 9,. 30952796,. 34 95 48 64 , 

.3 91 6 48 66,. 43533264,. 48029511,. 5263 5 06 9,. 57339 78 9,. 62 14 0105, 

.169 3 6998,. 1693 6998,. 185645 67,. 214811 81,. 25 04 91 6,. 28 990 367,. 331 73 77 5, 

.375 27812,. 4 201 0459,. 4 65 96 496,. 5127127 7,. 560 27 373,. 60862959, 

. 12931031,. 129310 31,. 15261801,. 18711 985,. 22663938,. 26 894684,. 31300625, 
.35827161, .404 43 92 5,. 45133986,. 49 88 85 23,. 5 47 04 135,. 595 81 671, 
.10660 326,. 106603 26,. 13 49 8 60 8,. 17 27 9866,. 2 14 6 26 72,. 25860856,. 30 390795, 
.350 10 593,. 39 69 8293,. 44 442 55 5,. 49 238 31 8,. 5 40 84699,. 5898 4179, 
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DLSOR Cl, 1 ) = 0. „ 0., 0 . , .0000002, . 000 005 7 3, . 0000 52 36 , . 0 00 26419, . 000 901 25 , 
.0023587, .00514317/ .00989938, .00989938, .00989938, 

0.,0.,0., .0000031, .00005746, .00040704, .00162363, .00452431, .00996314, 
.01866783, .03116991, .04781193, .06880882, 

0., 0., .0000002 8, .0000 2575, .0 0035 647, . 001958 45, .006 40 297, 

.0153088 9, .02991039, . 05 091317, .0 78532 34, .11260453, .152 71789, 


O.,0.,. 00000 21 8,. 000 13 394,. 00140073,. 00623 07 6, .0174901, .0373587, 
.06695116,. 10642165,. 15522653,. 21239 05 4,. 2 767 23 71, 

0., 0.,. 00001098,. 00048559,. 00397 466,. 01496156,. 0373917 2,. 0 7337 245, 

.12323 764,. 186 0700 2,. 26023 292,. 3437 6991 ,. 4346 7831, 

0.,0.,. 00 00 4119,. 00 134715,. 00899 29 ,.0 29 6 67 45,. 06773287, 

.1244 23 98,. 19867459,. 28812978,. 38 990179,. 50105232,. 61 88581, 

0.,. 00000012,. 0001 2344,. 003 0712,. 01732 98 3,. 05143 03,. 10922 78 7, 

.19030 35 9,. 2 91 7307 9,. 4095372,. 5 395 9966,. 67812008,. 82183021, 

0.,. 00 00002,. 0003097,. 006 05 144,. 02 97 10 67,. 0808637 4,. 16185 55 50,. 26996484, 
.4 00 1399, .5467 85 65, .70468605, .86938119, 1.03 72605, 

0.,. 00000 18 5,. 001 31 425,. 017 3 06 57,. 06852806,. 16323 28 6,. 29 803 005, 

.46440 851, .65 298 40 8, . 8 5541135,1.0649734,1.2765648,1 .486 461, 

0., .00003844, .00889343, .06884447, .20612165, .41208982,. 66585564, 
.94742268,1.2414942,1.5375068,1.8285879,2.1104825,2.3807193, 

0.,. 000300 54,. 029811 73,. 16533512,. 41495684,. 74144 43 5, 1.107 25 41, 1.4860364, 
1.8616789,2.2251066,2.5716249,2.8991281,3.2070205, 

MT2(1) = 0.,. 2, .4, .6, .8,1. ,1.2, 1.4, 1.6, 1.8, 2. ,2. 2, 2. 4, 2. 6, 2. 8, 3. ,3. 2, 3. 4, 3. 6, 

3. 8. 4. . 4. 2. 4. 4. 4. 6. 4. 8. 5. .5. 2. 5. 4. 5. 6. 5. 8. 6. .6. 2. 6. 4. 6. 6. 6. 8. 7. .7. 2. 7. 4, 

7. 6. 7. 8. 8. . 8. 2. 8. 4. 8. 6. 8. 8. 9. .9. 2. 9. 4. 9. 6. 9. 8. 10., 

0=7., 

THODLT ( 1 ) = 1 .,.09695, .09616, . 09 487, .09315, .091 04, .08862, .08595, 

.08311,. 08 01 6,. 07714,. 0741,. 07 10 8,. 06812,. 065 22,. 06 24, 


. 05969, . 05708, 
.04374, .04189, 
.03262, .03135, 
.02406, .02321, 
. 01831, . 01773, 
I TER=2, 


05458. . 05219. . 04992. . 04775. . 0457, 

04013. . 03846. . 03688. .03539. .03397, 

03014. . 02899. .0279. .02686. .02588. .02495, 

02241.. 02164.. 02 091.. 020 21.. 01955.. 01 891, 

01717. . 01665. . 01614, 


T0 = 5 70 . , 
N = 5 2, 


X<1)«. 25, .5,1. ,2. ,3. ,4. ,5. ,6. ,7. ,8. ,9. ,10. ,11. ,12. ,13. ,14. ,15. ,16. ,17. ,18. ,19., 
20 . , 21 . ,22. , 23. , 24 ., 25 . ,26., 27 ., 28 . ,29. ,30. , 31 . ,32. ,33. ,34. ,35. ,36. ,37. ,38., 

39. ,40. ,41. ,42. ,43. ,44. ,45. , 46. ,47. ,48. ,49. ,50., 

Y<1)=. 0331, .066, .130 8, .2573, .3795, .4973, .6107, .7199, .82 47, .925 1,1. 0212, 1.113, 
1.2004,1.2835,1.3623,1.4367,1.5067,1.5725,1.6339,1.6909,1.7436,1.792,1.836, 
1.8757, 1.9111, 1.9421, 1.9687, 1.9911, 2. 0091, 2. 0227, 2. 032, 2. 037, 2. 0376, 2. 0339, 
2.0259,2.0135,1.9967,1.9757,1.9503,1.9205,1.8864,1.848,1.8052,1.7581,1.7067, 
1.6509,1 .5907, 1.5 263, 1.45- 75; 1.58 43, 1.3068, 1.225, 

RN=4 . 17E06, 


AM =2 . 2, 
$ 
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It should be noted that the input tables of conical flow parameters were obtained 
from reference 7 and need not be changed from case to case unless better coverage is 
desired. Similarly the input tables of 9/6 were obtained from reference 13 for a 
velocity -profile -index value of 7 and need not be changed unless better coverage is desired 
or the user wishes to examine the results for a value of the velocity profile index other 
than 7. Only odd values of the velocity profile index can be used in the program. 

Description of Output Data 

The output format of the program is arranged so that each iteration of a single case 
contains a listing of body coordinates, free-stream conditions, starting cone conditions, 
longitudinal variation of the inviscid parameters, shock coordinates, characteristic 
angles, longitudinal variations of the viscid parameters including flat-plate friction 
coefficients, integrated total -friction coefficients (body of revolution and flat-plate values) 
based on body surface area, and wave-drag, friction-drag, and total-drag coefficients 
based on body maximum cross-sectional area. It should be noted that the body coordi- 
nates presented for the second and succeeding iterations are the coordinates of the 
effective body and are used only in the solution of the inviscid parameters. The units 
shown in the output listing are correct only when the body coordinates are in inches. 

Example Output Listing 

The following is an example output listing for the RM-10 body at the conditions 
specified by the example input case previously presented: 
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RM 1 


0 BODY 



INPUT DATA 

ITERATION 


M= 2.20 



RN= 4.17000000E+06 

X 

Y 

DELTA 

.2500 

.0331 

7.4219 

1 .0000 

.1308 

7.4219 

2. Cuou 

.2573 

7. 2096 

3.0000 

.3795 

6.9670 

4.0UOG 

.4973 

6.7185 

5 . COuO 

.6107 

6.4697 

6 . OouG 

.7199 

6. 2320 

7.0000 

.8247 

5.9828 

8 . G U 0 0 

.9251 

5. 7333 

9.0000 

1.0212 

5.4893 

l o . Co o c 

1.1130 

5.2451 

11.0000 

1.2004 

4. 9950 

12. CuOO 

1.2835 

4.7504 

13.0000 

1.3623 

4. 5056 

14.0000 

1.4367 

4. 2550 

15.0000 

1.5067 

4.0042 

16 . 000o 

1.5725 

3.7646 

17.0000 

1.6339 

3. 5135 

18.000c: 

1 .6909 

3.2623 

1 9 . 0000 

1.7436 

3.0167 

20 . CGCC 

1.7920 

2.7710 

21 .0000 

1.8360 

2.5194 

22 . OOuO 

1.8757 

2.2734 

23.0000 

1.9111 

2.0274 

24.0000 

1.9421 

1. 7756 

25.0000 

1.9687 

1.5237 

26.0000 

1.9911 

1. 2832 

27.0000 

2.0091 

1.0312 

28. COCO 

2.0227 

. 7792 

29. OoOO 

2.0320 

. 532 8 

30 . 0000 

2 .0370 

. 2865 

31 . 0000 

2.0376 

.0344 

32.0000 

2.0339 

-.2120 

33.0000 

2.0259 

-.4584 

34 . 0000 

2.0135 

-.7104 

35 • 0000 

1. 9967 

-.9625 

36.0000 

l. 9757 

-1.2030 

37 • 0000 

1.9503 

-1.4550 

38.0000 

1.9205 

-1.7069 

39. 0000 

1.8864 

-1. 9530 

40 . OoOO 

1 . 8480 

-2.1991 

41. 0000 

1.8052 

-2. 4508 

42 . OOOO 

1.7581 

-2.6966 

43 . OoOO 

1 . 70 67 

-2. 9424 

44 . OOOO 

1.6509 

-3.1938 

45 .0 000 

1 .5907 

-3.4450 

46 • OOuO 

1 . 5263 

-3.6848 

47.0000 

1.4575 

-3.9357 

48. OOuO 

1.3843 

-4. 1866 

49. OoOO 

1.3068 

-4.4316 

50.0000 

1.2250 

-4. 6764 

STARTING CONE - 

2.07766 THET A= 

27.963 DELSO 


.0000 
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APPENDIX B 


Continued 


OUTPUT 


INV1SCI0 SOLUTION 


X ♦ I N 

P/P< INF) 

M 

CP 

X- SHOCK 

Y-SHOCK 

CHARACTERISTIC ANGLE 

D/Q,SQ IN 

.1250 

1.21015 

2.0777 

.0620 

.7569 

.4040 

36.1930 

.000213 

1ST OROER SHOCK 
.6250 

EXPANSION 

1.21015 

2.0777 

.0620 

3.50 98 

1.9672 

36.1930 

.003334 

1.5000 

1 .20054 

2 » j 828 

.0592 

7.3956 

4.1636 

35.9036 

.012463 

2. 5000 

1.18891 

2.089U 

.0558 

11.5841 

6.5177 

35.5675 

.026095 

3.5000 

1. 17743 

2.0952 

.0524 

16.1037 

9.0438 

35.2263 

.043088 

4 .5000 

1.16613 

2.1014 

. 0490 

20.9907 

11.7602 

34.8862 

.062444 

5.5000 

1.15557 

2 . 1 0 72 

.0459 

26.2378 

14.6616 

34 .5628 

.083405 

6.5000 

1. 14430 

2.1135 

.0426 

32.0112 

17. 8363 

34.2220 

. 105064 

7.5000 

1.1 3318 

2.1197 

.0393 

38.3270 

21.2900 

33.8818 

.126760 

8.5000 

1.12248 

2.1258 

.0362 

45. 2C99 

25.0338 

33.5503 

.148002 

9.5000 

1.11181 

2.1319 

.0330 

52.7965 

29.1381 

33.2186 

.168315 

10.5000 

1.10084 

2.1 382 

.0298 

61.2908 

33.7080 

32.8783 

.187222 

11 . 5000 

1.09039 

2.1443 

.0267 

70.6541 

38.7186 

32 . 5475 

. 204521 

12.5000 

1.08031 

2.1 503 

.0237 

81.1047 

44.2823 

32 .2194 

.220048 

13.5000 

1.07029 

2.1562 

.0207 

93.0007 

50.5827 

31.8856 

. 233621 

14.5000 

1 . 060o7 

2.1620 

.0179 

106.4822 

57.6874 

31.5548 

.245212 

15.5000 

1.05203 

2.16 72 

.0154 

1 21.4003 

65.5138 

31 .2432 

.254988 

16.5000 

1.04291 

2.1 728 

.0127 

139.1018 

74.7560 

30.9150 

.262822 

17.5000 

1.03401 

2.1783 

.0100 

159.728^ 

85.4751 

30.5899 

.268799 

18.5000 

1.02563 

2.1o35 

.0076 

1 83.5765 

97.8133 

30.2738 

.273101 

19.5000 

1.0173d 

2 . 1 8 8o 

.0051 

211.9466 

112.4267 

29.9584 

.275859 

20.5000 

1.00890 

2.1 940 

.0026 

246.9182 

130.3589 

29.6352 

.277177 

21 . 5000 

1.00092 

2.1991 

.0003 

289.0698 

151.8799 

29.3215 

.277302 

22.5000 

,9931o 

2.2040 

-.0020 

341 .6988 

178.6378 

29.0095 

.276452 

23.5000 

.98538 

2.2091 

-.0043 

410.5317 

213.4864 

28.6914 

.274833 

24.5000 

.97797 

2.2139 

-.0065 

501.5329 

259.3719 

28,3761 

.272707 

25.5000 

.97155 

2.2181 

-.0084 

619.9377 

318. 8651 

28.0805 

.270367 

26.5000 

.96475 

2.2 226 

- • 0 1C 4 

798 .0145 

408.0062 

27.7702 

.2680 14 

27.5000 

.95829 

2.2269 

-.0123 

1073.9261 

545.6279 

27.4625 

.265893 

28.5000 

.95244 

2.2308 

-.0140 

1534.3313 

774.5263 

27.1656 

.264230 

29.5000 

.94683 

2 .2 346 

-.0157 

2460.9632 

1233.7786 

26.8708 

.263227 

30.5000 

.94113 

2 . 2 3 84 

-.0174 

5189. 135C 

2581.5988 

26.5694 

.263093 

31 . 5000 

• 9 3 6l* 5 

2.2419 

-.0189 

34067.4106 

16807.7049 

26.2788 

.263987 

32.5000 

.93 1 2 1 

2.^452 

-.0203 

*2658.9214 

*1060.8293 

25.9903 

.266058 

33.5000 

. 92636 

2.2485 

-.0217 

-6264.6954 

-3028.7916 

25.6959 

.269479 

34.5000 

.92185 

2.2517 

-.0231 

-3437. 1813 

-1647.0495 

25.4045 

.274361 

35.5000 

.91830 

2.2541 

-.0241 

-2378.9731 

-1130.9714 

25.1329 

.28068 1 

36.5000 

. 9 1 444 

2.2568 

-.0253 

-1777. 1695 

-838.1242 

24. 8471 

.288593 

37.5000 

.91093 

2.2593 

-.0263 

-1414.4816 

-661.9880 

24.5645 

.298120 

38.5000 

• 90 809 

2.2613 

-.0271 

-1181.3130 

-548.9363 

24 . 2934 

.309183 

39.5000 

.90558 

2.2 630 

-.0279 

-1015.2246 

-468.5248 

24.0252 

.321738 

40.5000 

.90309 

2.2648 

-.0286 

-887.4902 

-406.7710 

23. 7516 

.335789 

41 . 5OC0 

.90125 

2.2 661 

-.0291 

-792.8251 

-361.0532 

23.4895 

.351157 

42.5000 

.89975 

2.2672 

-.0296 

-717.8015 

-324.8526 

23.2305 

.367713 

43.5000 

.89829 

2.2682 

-.0300 

-654.9294 

-294.5410 

22.9662 

.385383 

44.5000 

.89721 

2.2690 

-.0303 

-603.4036 

-269.7151 

22.7055 

.403983 

45.5000 

.89710 

2.2 690 

-.0304 

-563.6214 

-250.5486 

22.4648 

.423136 

46.5000 

.89674 

2.2693 

-.0305 

-527.0375 

-232.9267 

22.2106 

.442793 

47.5000 

.89680 

2.2693 

-.0305 

-495.8912 

-217.9235 

21.9603 

.462700 

48.5000 

. 8976U 

2.2687 

-.0302 

-470.2125 

-205.5484 

21.7224 

.482503 

49.5000 

. 89883 

2.2 6 78 

-.0299 

-447.8936 

-194.7846 

21.4885 

.501931 
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APPENDIX B - Continued 


VISCOUS SOLUTION 


X/L 

THETA, FT 

800Y-CF, LOCAL 

FLAT PLATE CF, LOCAL 

DEL STAR, FT 

DELTA, FT 

.00250 

.00001322 

.00467920 

.00347823 

.00004234 

.000 16715 

.01250 

.00005205 

.00396972 

.00322309 

.00016679 

.00066351 

.03000 

.00011052 

.00348725 

.00283543 

.00035526 

.00141481 

.05000 

.00017190 

.00322033 

.00262555 

.00055459 

.00220648 

.07000 

.00023040 

.00304933 

.00249651 

.00074600 

.00296303 

.09000 

.00028712 

.00292387 

. 00240630 

.00093298 

.00369826 

.11000 

.00034262 

.00282513 

.00233726 

.00111705 

.00441882 

. 1 300C 

.00039729 

.00274246 

.00228158 

.00129995 

.00513029 

. 15000 

.00045144 

.00267086 

.00223510 

.00148240 

.00583615 

. 1700C 

.00050524 

.00260815 

.00219531 

.00166484 

.00653841 

. 19000 

. U00558B5 

.00255243 

.00216059 

.00184790 

.00723914 

. 2 1C0C 

.0006 1244 

. 00250111 

.00212986 

.00203239 

.00794076 

.23000 

. 00066612 

.00245368 

.00210234 

.00221818 

.00864398 

.25000 

, OuO 7 1997 

.00241038 

.00207743 

.00240558 

.00934994 

.27000 

.00077415 

.00236987 

.00205471 

.00259533 

.01006067 

.29000 

.00082877 

.00233093 

.00203387 

.00278750 

.01077719 

. 31000 

.00088387 

.00229607 

.00201460 

. GO 2981 57 

.01149899 

.33000 

.00093959 

.00226351 

.00199671 

.00317944 

.01222969 

.35000 

.00099607 

.00223117 

.00198003 

.00338091 

.01296997 

.37000 

.00105337 

.00220080 

. 00196441 

.00358581 

.01372001 

. 39000 

.0011 1159 

.00217295 

.00194973 

.00379492 

.01448152 

.41000 

• u0117o9J 

.00214525 

.00193588 

.00400931 

.01525703 

.43000 

.00123138 

.00211806 

. 00192279 

.00422831 

.01604618 

. 45000 

.00129315 

.00209289 

.00191038 

.004452 65 

.01685062 

.47000 

.00135641 

.00206817 

.00189858 

.00468343 

.01767307 

.49000 

.00142132 

• C0204205 

.00183735 

.00492060 

.01851447 

.51000 

.00148792 

.00201914 

.00187663 

.00516306 

.01937344 

.53000 

.00155653 

.00199795 

.00186638 

.00541436 

.0202566 2 

.55000 

.00162737 

.00197497 

.00185657 

.00567397 

.021 16451 

. 57000 

.00170059 

.00195311 

.00184715 

.00594179 

.02209778 

. 59000 

.00177646 

.00193372 

.00183810 

.00621940 

.02305979 

.61000 

.60185334 

.00191353 

.00182940 

.00650889 

.02405529 

.63000 

.00193748 

.00189306 

.00182102 

.00680949 

.02508433 

.65000 

.00202325 

.00187500 

.00181294 

.00712320 

.02615116 

.67000 

.U0211315 

.00185672 

.00180513 

.00745260 

.027 261 6 5 

. 69000 

.00220768 

.00183527 

.00179759 

.00779833 

.02841860 

.71000 

. 00230718 

.00181726 

.00179029 

.00815962 

.02962197 

, 73000 

.00241250 

.00180126 

.00178322 

.00854331 

.03088426 

. 75000 

.00252442 

.00178152 

.00177637 

.00894998 

.032 20 894 

. 77000 

.00264369 

.00176237 

.00176973 

.00938103 

.03359973 

. 7900C 

.00277136 

.00174633 

.00176328 

.009841 16 

.03506568 

.81000 

.00290881 

.00172848 

.00175701 

.01033631 

.03661812 

. 8300 C 

.00305740 

.00170953 

.00175091 

.01086853 

.03826203 

. 85000 

.00321893 

.00169396 

. 00174498 

.01144524 

.04000961 

. 8700 0 

.60339580 

.00167812 

.00173921 

.01207592 

.041 87622 

, 89000 

.00359079 

.00165649 

.00173358 

.01276855 

.04387364 

.91000 

. 00380694 

.00164002 

.00172810 

.01353017 

.04600892 

, 93000 

.00404904 

.00162747 

.00172275 

.01438348 

.04830911 

,95000 

.00432277 

.00160874 

.00171753 

.0 1534420 

.05078633 

, 97000 

.00463536 

.001591 96 

.00171243 

.01643445 

.05345139 

,99000 

.00499670 

.00106411 

.00170962 

.01768930 

.05631745 


BODY— CF»T3TAL= .00200549 FLAT-PLATE CF,TOTAL= .00189077 

CD,WAVE= .03848177 CD,FRICT10N= .07321805 CD *TOT AL= .11169982 
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APPENDIX B - Continued 


RM 10 BODY 


I N PUT DATA ITERATION 

M= 2.20 

RN = 4. 17000000E+06 


X 

Y 

DELTA 

.2500 

.0 340 

7.5772 

1.0000 

. 1337 

7.5772 

2.0000 

.2627 

7. 3486 

3.0000 

. 3872 

7. 0981 

4. 0000 

.5073 

6.8455 

5.0000 

.6229 

6. 5944 

6.0000 

.7343 

6.3554 

7.0000 

. 8412 

6. 1057 

8.0000 

.9438 

5.8561 

9.0000 

1.0421 

5. 6124 

10.0O0G 

1 .1360 

5.3689 

11.0000 

1.2256 

5. 1197 

12.0000 

1 .3109 

4. 8762 

13.0000 

1 . 3920 

4.6327 

14. 0000 

1.4686 

4. 3837 

15.0000 

1 . 5409 

4.1344 

16.0000 

1.6090 

3. 8967 

17.0000 

1.6728 

3.6481 

18. 0000 

1.7322 

3. 3992 

19.0000 

1.7873 

3. 1561 

20.0000 

1 . 8382 

2. 9135 

21.0000 

1.8848 

2.6651 

22.0000 

1.9271 

2.4224 

23.0000 

1.9652 

2.1802 

24. 0000 

1 .9989 

1. 9325 

25.0000 

2.0233 

1.6844 

26 . OOOO 

2.0536 

1.4484 

27.0000 

2.0746 

1.2019 

28.0000 

2.0912 

.9550 

29.0000 

2.1037 

. 7144 

30.0000 

2.1120 

.4748 

31.0000 

2. 1160 

. 2299 

32.0000 

2.1158 

-.0090 

33.0000 

2.1115 

-. 2464 

34.0000 

2 .1030 

-.4887 

35 . 0000 

2.0903 

-. 7310 

36.0000 

2.0735 

-.9601 

37.0000 

2.0526 

-1.1984 

38 . 0000 

2.0275 

-1.4363 

39.0000 

1.9984 

-1. 6667 

40.0000 

1.9653 

-1.8940 

41.0000 

1.9282 

-2. 1249 

42.0000 

1.8872 

-2. 3476 

43.0000 

1.8424 

-2.5657 

44.0000 

1.7938 

-2. 7850 

45 . OOOo 

1.7414 

-3.0010 

46. 0000 

1.6855 

-3. 1982 

47.0000 

1.6261 

-3.3976 

48.0000 

1 . 5634 

-3. 5892 

49.0000 

1 .4976 

-3.7633 

50.0000 

1.4280 

-3.9830 


STARTING CONE - M= 2.07391 THET A= 28.003 DELSOR= 


.0000 
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APPENDIX B 


Continued 


OUTPUT 

INVISCID SOLUTION 


X, IN 

P/P( INF) 

M 

CP 

X- SHOCK 

Y-SHOCK 

CHARACTERISTIC ANGLE 

D/0, SO IN 

. 1250 

1.21727 

2.0739 

.0641 

.7452 

.3991 

36.4C51 

.000221 

1ST ORDER SHOCK 

EXPANSION 







.6250 

1 . 2 1 72 7 

2.0739 

.0641 

3.4495 

1.9400 

36.4C51 

.003447 

1 . 5000 

1 . 20 642 

2.0797 

.0609 

7.2787 

4.1105 

36.0895 

.012844 

2.5000 

Z. 19412 

2.0862 

.0573 

11.4037 

6.4349 

35.7402 

.026851 

3.5000 

1.18245 

2.092b 

. 0539 

15.8461 

8.9239 

35.3937 

.044326 

4.5000 

1.17105 

2.0987 

.0505 

20.6422 

11.5962 

35.0504 

.064254 

5.5000 

1. 16042 

2.1 045 

.0473 

25.7831 

14.4457 

34.7254 

.085869 

6.5003 

1 . 14912 

2.1108 

.0440 

31.4263 

17.5562 

34.3839 

.108252 

7.5000 

1.13798 

2 . 1 1 70 

.0407 

37.5866 

20.9330 

34.0436 

. 130729 

8.5003 

1.12727 

2.1231 

.0376 

44.2335 

24.5843 

33.7124 

.152803 

9. 5C00 

1.11663 

2.1292 

.0344 

51.6422 

28.5750 

33.3818 

. 173992 

10 .5000 

1.10568 

2.1355 

.0312 

59.8560 

33.0049 

33.0427 

.193806 

11.5000 

1.09522 

2.1415 

.0281 

68.8852 

37.8487 

32.7129 

.212030 

12.5000 

1 .08499 

2 . 1 4 75 

.0251 

78.9395 

43.2143 

32.3850 

.228461 

13.5000 

1.07480 

2.1536 

. 0221 

90. 3357 

49.2641 

32.0516 

.242905 

14.5000 

1.C6501 

2.1 594 

.0192 

103.2005 

56.0591 

31 .7209 

.255325 

15.5000 

1.05626 

2.1647 

.0166 

117.3527 

63.5000 

31.4103 

.265895 

16.5000 

1.047J7 

2.1 703 

.0139 

134.0140 

72.2182 

31 .0849 

.274488 

17.5000 

1.03812 

2.1758 

.0113 

153.2984 

82.261 5 

30.7610 

.281186 

18.5000 

1.02972 

2.1810 

.0088 

175.4076 

93.7250 

30 .4473 

.286174 

19.5C00 

1.02149 

2 , 1 8ol 

.0063 

201.4188 

107.1527 

30.1353 

.289583 

20.5000 

1.01303 

2.1914 

.0038 

233.1194 

123.4436 

29.8156 

.291513 

21.5000 

1.005U6 

2.1964 

.0015 

270.8507 

142.7509 

29.5054 

.292204 

22.5000 

.99728 

2.2014 

-.0008 

317.2018 

166.3695 

29.1971 

.291865 

23.5000 

.98944 

2.2065 

-.0031 

376.6586 

196.5378 

28.8827 

.290696 

24.5000 

.98191 

2.2113 

-.0053 

453.5457 

235.3904 

28.5702 

.288951 

25.5003 

.97539 

2.2156 

-.0073 

550.3505 

284.1334 

28.2785 

.286927 

26.5000 

.96855 

2.2 201 

-.0093 

689.4691 

353.9190 

27.9732 

.284827 

27.5000 

• 9 62C i 

2.2 244 

-.0112 

892.7817 

455.5385 

27.6702 

.282895 

28.5000 

.95610 

2.2 2 8*+ 

-.0130 

1202.6105 

609.8897 

27.3786 

.281360 

29.5000 

.95049 

2.2321 

-.0146 

1737.4932 

875.5278 

27.0906 

.280426 

30.5000 

.94462 

2.2359 

-.0163 

2889.6385 

1445.9035 

26.7966 

.280301 

31.5000 

.93972 

2.2394 

-.0178 

6363.0721 

3160.5373 

26.5135 

.281143 

32.5003 

.93490 

2.2427 

-.0192 

84627.1568 

41690.3154 

26.2342 

.283104 

33 .5000 

> 9 30U6 

2.2 4o0 

-.0206 

*6249 ,5759 

*7655.2253 

25.9498 

.286352 

34.5000 

.92551 

2.2491 

-.0220 

*2140. 1962 

-5848.9846 

25.6677 

.291006 

35.5000 

.92197 

2.2516 

-.0230 

-7228.2274 

-3448.4424 

25.4C78 

.297042 

36.5000 

.91819 

2.2542 

-.0241 

-4974.0266 

-2349.1607 

25.1363 

.304607 

37.5000 

.91472 

2.2566 

-.0252 

-3763.5372 

-1760 .3202 

24.8682 

.313728 

38.5003 

.91197 

2.258o 

0260 

-3044. 5403 

-1410.6392 

24.6136 

.324325 

39.5000 

. 90962 

2,2602 

-.0267 

-2561 .9867 

-1 1 76.4680 

24.3657 

. 336343 

40.5003 

.90734 

2.2618 

-.0273 

-2204. 7136 

-1003.3193 

24.1147 

.349777 

41 . 5000 

.90576 

2.2629 

-.0278 

-1948.7660 

-879.3893 

23.8781 

.364443 

42 .5000 

• 90464 

2.2637 

-.0281 

-1753.5097 

-784.9071 

23.6502 

. 380190 

43.5000 

.90367 

2.2644 

-.0284 

-1594.4185 

-707.9681 

23.4222 

.396926 

44.5000 

.90309 

2.2648 

-.0286 

-1466.4897 

-646.1202 

23.2012 

.414462 

45.5000 

.90372 

2.2644 

-.0284 

-1373.3165 

-601.0587 

23.0096 

.432383 

46.5000 

.90433 

2. 2o39 

-.0282 

-1291.1833 

-561.3225 

22.8156 

.450593 

47. 5000 

.90550 

2.2631 

-.0279 

-1224.1866 

-528.8876 

22.6342 

.468822 

48.5000 

. 90768 

2.2616 

-.0273 

-11 74.0993 

-504.6086 

22.4793 

.486677 

49.5000 

.90751 

2.2617 

-.0273 

-1 1C7.2445 

-472.2049 

22.2581 

.504437 
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APPENDIX B - Concluded 


VISCOUS SOLUTION 


X/L 

THETA, FT 800Y-CF , LOCAL FLAT 

PLATE CF, LOCAL 

DEL STAR, FT 

DELTA, FT 

00250 

• 00 00 1323 

.00469244 

.00347699 

.00004228 

.00016717 

01250 

.OGUO5210 

.00398097 

.00322194 

.00016659 

.00066362 

03000 

.00011062 

.00349524 

.00283449 

.00035489 

.00141508 

05000 

.UU017203 

.00322666 

.00262477 

.00055405 

.00220683 

07000 

•Gu023057 

.00305475 

.00249583 

.00074525 

.00296335 

09000 

.00028732 

.00292901 

.00240568 

.00093198 

.00369851 

1 1000 

.00034284 

. C0282925 

.00233669 

.00 111580 

.00441898 

1 3000 

. GOU39754 

.00274658 

.00228105 

.00129838 

.00513025 

15000 

.00045171 

.00267488 

.00223461 

.00148052 

.00583590 

1 700 C 

.00050554 

.00261171 

.00219485 

.00166262 

.00653792 

19000 

.00055918 

.00255549 

.00216015 

.00184529 

.00723832 

21000 

.00061278 

.00250456 

.00212944 

.00202933 

.00793950 

23000 

. 00066o48 

.00245701 

.00210194 

.00221473 

.00864238 

25000 

.00072036 

.00241337 

.00207705 

.00240183 

.00934811 

2700 0 

.00077456 

.00237289 

.00205436 

.00259123 

.01305854 

2900^ 

.00082920 

.00233439 

.00203353 

.00278313 

.01077484 

31000 

.00088432 

.00229778 

.00201429 

.00297701 

.01149656 

33000 

. 0 0o94o04 

.00226543 

.00199641 

.00317431 

.01222657 

35000 

.00099653 

.00223335 

.00197975 

.00337533 

.01296633 

3700 r 

.00105384 

.00220239 

.00196415 

.00357981 

.01371594 

39C00 

.00111207 

.00217380 

.00194948 

.00378832 

.01447671 

41000 

.00117137 

.00214683 

.00193565 

.004001 90 

.01525116 

43000 

.00123186 

.00211946 

.00192258 

.00422037 

.01603968 

45000 

.00 129364 

.00209389 

.00191018 

.00444409 

.01684336 

4700C 

.00135689 

.00206967 

.00189841 

.00467415 

.01766490 

49000 

.00142181 

.00204485 

.00188720 

.00491089 

.01850585 

51000 

.00148845 

.00201986 

.00187649 

.00515342 

.0193651 8 

53000 

.00155704 

.00199903 

.00186626 

.00540372 

.02024699 

55 00 0 

.00 162788 

.00197690 

.00185647 

.00566284 

.021 15436 

57000 

.00170113 

.00195451 

.00184708 

.00593047 

.02208757 

59000 

.00177700 

.00193403 

.00183805 

.00620746 

.02304887 

61000 

.00185586 

.00191494 

.00182937 

.00649584 

.02404287 

6300 C 

,001938o2 

.00189437 

.00182102 

.00679612 

.02507172 

65000 

.00202379 

. C01 87 542 

.00181296 

.00710910 

.02613772 

67000 

.00211367 

.00185783 

.00180519 

.00743728 

.02724663 

69000 

.062^0821 

.00183878 

.00179768 

.00778245 

.02340309 

71000 

. U02307 78 

.00181799 

.00179041 

.00814440 

.02960794 

73000 

.00241306 

.00180221 

.00178339 

.00852609 

.03386755 

75000 

.00252499 

.00178455 

.00177658 

.00893202 

.032 19158 

77000 

» 00 264431 

.00176500 

.00176998 

.00936313 

.03358305 

79000 

.00277200 

. 00174720 

.00176358 

.00982232 

.03504818 

81000 

.00290941 

.00173115 

.00175737 

.01031527 

.03559799 

83000 

.00305805 

.00171225 

.00175134 

.01084704 

.03824205 

85000 

.00321958 

.00169442 

.00174549 

.01142187 

.03998779 

87000 

.00339636 

.00167904 

.00173981 

.01204882 

.04185304 

89000 

.00359131 

.00166161 

.00173428 

.01273842 

.04384447 

91000 

.00380756 

.00163933 

.00172891 

.01349937 

.04598047 

93000 

.00404939 

.00162515 

.00172371 

.01434410 

.04827087 

95000 

.00432296 

.00160982 

.00171865 

.01529773 

.05374130 

97000 

.00463544 

.00159164 

.00171375 

.01638138 

.05340133 

99000 

.00499657 

.00107464 

.00171112 

.01762698 

.05626034 


BODY-CF, TOT AL= 

.00200748 

FLAT-PLATE CF, 

TOTAL= .00189082 

CD, 

WAVE= .03867395 

CD, FRICTION 

= .07329059 

CD, TOT AL = 

.11196454 
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Figure 5.- Comparison of the theoretical methods of reference 5 (eq. (2)) with the present theory (eq. (3)) 
for estimating the momentum thickness of a turbulent boundary layer. 







Second iteration 
First iteration 




(a) Inviscid solution. 


Figure 8.- Typical results of the present method applied to an Adam's body of FR = 10 and 

at Mo, = 4.0 and R^j = 9.0 x #. 
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Figure 10.- Variation of momentum thickness and local skin-friction coefficient on NACA RM-10 missile body. 
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Figure 11.- Comparison of theoretical and experimental friction-drag and total-drag coefficients 
for NACA RM-10 missile body. 







Figure 13.- Comparison of theoretical and experimental friction-drag and total-drag coefficients 
for the waisted body. R ro _i = 10 x 10°. 
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Figure 14.- Evaluation of theoretical total-drag predictions. 
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